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A  hstmrt 

l  lir  Ix'liHvior  of  lliiids  iia-H  Ix-tMi  stiKlii**!  for  wars.  I'liis  Ix-liavior  lia,s  Ixs-n 
MKxIrIril  liy  lli«>  Nav  irr-Slok«‘s  r<|iialioiis  siixr  I  ho  ini^l-iiiiii'lis'iil  li  lio\V('\«'r. 

tlioso  )<«|iial  ions  arc  so  (oiiiploN  dial  (lioir  iiso  hy  onj^iixs'is  was  iiii|)ra«  I  i<  al  until  I  hr 
a<l\  rnl  <»r  I  hr  nitxlrrn  (  oinpiitri  .  I  hr  saiiir  i  oinpiilrrs  lhal  matir  I  lirsr  r(|nalions 
nsrfnl  rrralnl  a  nrw  prohlriii;  <lala  saliiralion.  I  lx*  sohilion  to  |Ik‘  Navirr-.Slok<'s 
rtpiHlions  which  rrpr«'s<Mil  an  at'rrxlyiiaitiic  flow  lirhl  is  cuinprisrd  of  Ihoiisainis  if 
not  millions  of  nnmix'rs.  This  mncli  data  makes  it  virtually  im|K>ssil)le  to 
coiK'rpI  nali/r  the  sit  nation  within  tin*  flow  field.  To  alleviate  this  prohlem,  various 
computer  graphics  teclinirpies  have  her'ii  used  to  provide  images  of  flow  fields. 

Some  of  tiu'se  |.<>chni<pi<‘s  r«'ly  on  approximating  isometric  snrfa«  «*s  of  interest  with 
polygons:  others  only  proviih*  an  imag<*  of  a  two-dimensional  “slice”  of  the  flow 
Ih'hl. 

I  his  thesis  projri  I  applies  a  relatively  new  graphics  ledinitpie  known  as 
“volumetric  r«‘n<lering”  to  generate  three-dimensional  imagr's  of  aiuodynamic  flow 
fields,  \olnmetric  rendering  has  shown  promising  rr'snlts  in  several  applications  of 
s<  i<Mitific  visualisation.  One  advantage  of  this  t<*chni(|ue  is  that  it  doesn't  rely  on 
gi'ometric  primitives  to  approximate  a  surface.  .Also,  no  n  priori  knowle»|gi'  of  the 
flow  held  s  form  is  rispiired  to  pro<lu<e  tiu*  final  image. 

To  assess  the  utility  of  volimwtric  rendering  in  the  field  of  computational 
fluid  d  ynamics  (Cl'l)).  a  program  known  as  \  lIM.H  is  developed.  This  thesis 
outlines  the  r<’<purem<'ut s  and  spe<  ifi»  algoiitluu''  hn  this  piogram.  Images  of  How 
fields  ate  presented  and  discuss<’d  ah>ug  with  pronram  I'lihancenient s.  .Mthoiigli 
\  Il’KIt  is  developed  foi  ('l  l)  purpose's,  its  desinn  permits  its  use  iii  dilfereiil 
applications,  sme  ral  of  whie  h  are  dise  uss*m|. 
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I.  IntrtMUiriion. 


I.l  Tht  •iis  Slalf  mmt 

A  r*'!aliv<'ly  now  oompuU'r  j^tapliics  tediiii«|iio  known  aa  volniiioliio 
|■on*l«M•ing  is  hoing  rfrognizo*!  as  a  poworful  means  of  sciontilir  visualization  in  sncIi 
ilivorso  fiokls  as  m<*l»sM'ology.  mo<lioino.  an«l  molonilar  pliysirs.  This  thosis  will 
show  t  hat  volumotri<-  roiuloring  shonhl  also  he  dovolopod  as  an  ongimsM'ing  tool  for 
pro<hiring  <l<'tail<'<l  images  of  th<>oroti<  ally-  ami  expoiimontally-rlorivod 
I hns’-dimonsional  ('U))  aorodynamic  (low  (ndds.  liv  providing  a  dl)  Image  of  this 
<lata.  the  time  r<’<piii(sl  to  t«*st  ami  |M-if<rl  mTodyuamio  designs  may  Im* 
signilioanlly  rednood.  thus  increasing  oiigimHuing  elfici«'nry  and  productivity. 

I.J  liackyround 

f  low  fiehl  v  isualization  is  an  ♦’ssential  tool  in  many  liclds  of  engimsuing. 
W’lu'ther  a  vlesign  is  for  a  huilding  that  must  withstand  high  winds,  a  ship  that 
must  negotiate  lu-avy  seas,  or  a  su|>ersonic  fighter  that  must  pass  through  the  air, 

I  he  eiigims'i  must  he  ahle  loslinlv  luiw  t he  st i u« )  me  iiilcrat  ts  with  its  fluid 
eiiv  iromm’iit .  f’.ver  since  engineers  starteil  workinu,  witli  lluid  dynamic  s  prohlc'uis. 
they  have  hc'en  faccsl  with  the  same-  prohlem:  how  to  \  isualize  the-  iuv  isilile.  i  util 
m)w.  c'liginec'rs  have  relic’cl  cm  two  l>asi<-  approac  lies  to  this  prohlem.  I  hc'  first  of 
tlu’se  approac  hes  is  to  huilci  a  phvsic  al  mocjc'l  ol  thc“  ohjc'cl  undcM  ch'sign  and 
actually  ohserve  how  the  ohjc'ct  behaves  in  a  flow  field.  This  technicpie  typically 
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ii|i|)r«>ii('li  is  U»  MSI'  H  ('oiii|)iitri  lo  si»|\<>  (III*  Navii*r-.S(ok<‘s  i'(|iial ions  for  (In*  olijrrt 
ami  ri'mliM'  a  I wo-ilininisional  (21))  iiiia|;<‘  liasisl  on  tliosc*  solnlions.  Moth  of  llii'si- 
a|>|>ioa»  lirs  liavr  si'vrn*  limitations. 

rin*  \N  I  ij;hl  In'otliors  wcri*  |iionis'is  in  tin*  ili'vi'lopini’nl  of  «’X|MTini«‘n(al  How 
lifhls.  I'lii'ii  lirsl  ilt'\i'lo|)m('nl  in  (liisaira  was  a  !>i<  \ (■|«‘  inonn(i‘il  ilov  ire  wliii  li 
lii'hl  moilrl  airfoils  in  llii*  How  lirhl  (hat  <lrv<*lo|H*<l  as  (hr  hiryi  Ir  was  prilalril 
Ihron^h  town.  .Alt hough  nsrfni.  this  arrangrmrnt  was  too  snsrrpt ihir  to 
niK'ont  rollril  (omiitions  (sm  h  as  variaUIr  win«ls)  ami  inaernratr  inrasnrrmriits  (fine 
to  operator  overload)  to  he  useil  for  tietailed  engineering  work.  To  overcome  these 
shortcomings  the  Wrights  invenleil  the  wind  tnimel.  a  device  that  conld  generate  a 
constant,  conirollefi  flow  fiehl  wherein  the  interaction  l)etwe<*n  airfoils  anil  their 
(Inid  environment  conltl  he  carefnlly  measnrefi  (  1:1-')  1).  Sinci*  then,  f'ligirn'ers  have 
relif’d  e.xtensively  on  flow  Held  stndif's  performed  in  wind  tnnni'ls.  Although  the 
wintl  tnmiel  lia.s  proven  itself  an  invalnahle  t<K)l  to  the  I'ngineer.  it  reipiires 
cnmhersome  ami  e.xpeiisive  |»rocednres  ( 1:1 ). 

do  begin  wit  h,  a  physical  mtxlel  of  the  ohject  under  test  (OUT)  ninst  he 
cf)iistrncted.  riiesf*  moflels  range  in  size  from  nmler  an  inch  in  length  to  full-size 
moilel  aircraft.  If  tin*  ohject  s  design  nei'fls  to  he  changeil.  the  mixlel  mnst  also  he 
he  MKxIilied  or  Cfxnpletf'ly  rehnilt.  When  the  mixh'l  is  rmished.  it  mnst  he  pri'pared 
for  tf’sting.  This  process  reipiires  monnting  (he  moilel  in  a  wind  tunnel  and 
preparing  (he  instrumentation  reipiired  for  the  test.  Once  (he  model  and  wind 
(nnnel  are  ready,  (he  test  is  run.  I'inallv.  after  the  test  is  completed,  the 
iiist  rnmeiit  at  ion  and  model  mnst  he  removed  Irom  the  tunnel.  II  more  data  is 
leipiired  later.  I  he  entire  test  seiinence  must  he  lepealed.  I'.very  step  ol  (his  process 
involves  time  ami  money;  in  an  era  of  tighteneil  resean  h  hndgets.  this  i  ycle  is  (ix) 
expensi\e  to  he  performed  often. 
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wiiiil  liniiK'Is.  I'linr  air  |■l•lHliv^ly  f«'\v  wiml  tiiiiiii-)s  availalili'  Id  niniiM'i'i'H.  If  (lino 
is  no  wiml  liiiinrl  availalilr,  an  niRiiusT  s  pi-Djis-l  niiisl  lx-  ilclayisl  or  possilily 
♦  aiu  i'ioil  allo^i'l  lior. 

I'orl  imali'ly,  lln'  ns  onl  iiicroaso  in  ruiiipiitor  (  apaliilitirs  aiul  availaliilily 
liavi'  pirsonlrd  solutions  lo  many  of  tlioso  prolilonis.  To  lio)>in  willi.  I  In*  compnli'i 
nsos  a  moilol  llial  is  alist  rarl  rallior  (lian  pliysical.  Wlini  rompanil  lo  physical 
nuxlols.  an  alistiart  moih'l  is  r<'la( iv<-ly  I'asy  lo  roi'sl rin  l  anil  modify.  Secondly, 
computers  are  now  commonplace  and  available  to  virtually  every  engineer.  Finally, 
no  speciali/eil  insirumentation  is  ni'tsieil  to  run  a  How  Pield  tesi  on  a  compnter; 
only  a  <  ompnt«'r  with  sulFicieiit  memory  and  some  means  of  displaying  the  U‘8t 
resnils  is  retpiiriMl.  llowi'ver,  lh<*  .sam<‘  computers  thal  provide  this  alwtract  design 
capabilily  also  creale  a  new  problem:  data  saturation.  'I’herefore,  computer 
availability  is  nol  enough;  a  satisfactory  method  for  displaying  ainl  analyzing  test 
r«’sults  still  poses  a  serious  problem  lo  if'searcliers. 

I'  low  held  data  is  normally  represeiiled  by  a  lhre<*-«limensional  (dl))  mesh  of 
data  points,  d  h«‘s«'  data  points  are  compose<l  of  solutions  to  Navier-Stokes 
erpiations  (a  set  of  mat luMiial ical  erpialions  that  model  fluid  behavior).  F.ach 
solution  may  contain  values  for  as  many  as  twenty  variables.  If  a  flow  field  is 
niodeh’d  by  a  lUU.x  lUO.x UK)  no<le  mesh  tln'ie  will  be  l.lHK).UtH)  locations  where  data 
is  to  b<-  analyz<*d.  This  volume  of  data  is  too  much  for  a  engimvr  to  assimilate 
from  printed  nnmIxMs  alone.  Therefore,  a  graphic  representation  is  necessary. 

F.arly  vdforts  at  flow  fudd  visualization  were  based  on  nuxleling  the  action  of 
wind  tunnels  themselves.  One  such  systv-m,  called  ('INFM.\.  was  ilesignevi  by  ('apt 
F.  P.  .Ambnrn  and  used  a  le«  hni<|n<‘  known  as  "parliv  le  lra<  ing".  This  system 
simnlalefi  the  release  of  parlich’s  in  a  lh»w  li«'ld  and  trati'd  thi'  part  ivies'  paths 
through  I Ix'  (iehl  as  a  function  ol  elapseil  lime  (see  l  ignre  1  ).  boeal  velocities  at 
each  node  wv're  stored  in  a  table.  As  a  particle  tjeared  a  node,  its  location  ainl 


\\<'ic  ii|ul.il<-i|  li\  ii'li'K'iH  iiii^  I  Ilf  iiimI)'  l.ilili'  I  I I  J'l).  ('INI'.MV  |M(mIii<  )'<l 
■Jl)  and  .'ll)  iiiiagrs  and  iis<*d  v«'«'(oi-lia.*i«i|  f^rapliits  syslcins  (I;l(i-I7).  lids  svsinn 
wa.s  an  ini|»i«»\«Mnrn(  «»vn  lonlriniioiarv  s\ sinus  in  llial  il  displayrd  l«'ni|><)ral 
fralnrrs,  ll«»wrvn\  tin*  nsr  *>f  a  svsinn  |»rn  lnilnl  a  anilinnoiis  display 

ImMwismi  data  puinis.  Tlins,  if  llir  nsrr  r«>l**asnl  paiiirlrs  in  (lie  wrtniK  areas,  rriliial 
fi'ainres  nmld  lie  oinillt'd  fnim  llie  linal  flow  field  iina^^e  {\nh:  \ ntor-liased 
displays  lia\«‘  tin'  advantage  of  <  lean.  «  risp  lines  and  i  iirves;  liowj'ver.  it  is  iliflii  till 
to  fill  in  solid  areas  (palelies)  of  ('olor  due  to  liardware  liming  constraints. 
Kaster-lia.st'd  graphics  displays  are  l>(‘lt<‘r  suitnl  for  palcli-Klling.  A  full  di.scus.sion 
of  vector-  and  raster-ha.s<'«l  graphics  displays  may  he  found  in  Chapter  .1  of  (7).).  In 
spile  of  this  shortcoming,  partich'  tracing  is  still  one  of  the  more  |>opular  technicpies 
iisnl  for  rendering  flow  fiehl  images.  Il  is  a  relatively  fast  l('chni(pie  and  lends  itself 
nicf'ly  to  user  interaction.  (IH APlI.'lI).  for  example,  is  a  software  package  develo|)etl 
for  N.'XSA's  Ann'S  |fesear<h  Cent<>r.  Although  si'veral  ilifferent  methods  for 
reiul<Miug  images  of  (low  (iehls  are  provuhxl.  the  particle  tracing  methorl  is  still 
consi<h'r«’<l  tin'  most  useful  te<  hui(|u<' of  those  available  on  (III  Al’ll-ll)  (2). 

In  ins.').  Ki’iineth  Kroos  pr«'sente<l  a  stinly  of  :ll)  visualization  techniipies 
involving  ve<  lor-hase<l  displays.  Kroos*  efforts  involveil  displaying  various  sha|)es  at 
<'a«  h  data  point  to  show  trends  within  a  flow  fiehl.  I'lie  primary  objective  of  tln’se 
•’xperimenis  wa.s  to  develop  an  algorithm  that  couhl  «piickly  general**  tln*se  images 
( 1  f:i;J()).  Kroos  concluih’<l  that  the  optimum  graphic  symlwl  for  this  purpose  was  a 
simph'  arrow-shaped  vector.  Such  a  vector  is  “am  hored"  to  each  rial  a  point  and 
has  its  magnitude  and  <lirection  d»’termim'<l  by  (In'  Iluid  flow  at  its  origin. 

.Mihough  this  lechni(|ue  allows  fast  remleriiiu,  of  an  iiuagt'.  Kroos  notes  that  il  has 
s*'\  <'ral  drawbacks  (11:111).  I'  irst .  du*'  lo  I  he  relal  ivel\  small  si/e  of  most  graphii  s 
(h'\  i<  *'s  il  is  «lilli*  lilt  lo  <lispla\  vi'iy  larg**  ami  \er\  small  \«'clois  in  I  In'  sam*'  image. 
Anolln'r  problem  pos«‘d  by  this  l<'<hni<(ue  is  its  inability  lo  ri'lh'i  l  li'inporal 
h'alures:  (hns  multiple  images  must  In'  pro«hn<'<l  to  show  t ime-cl('p<*mh'nt  features 
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Figure  1.  'rypiral  I’arlule  Tracing  Image 

williiii  I  lie  How  lii'UI.  Also,  when  reii<lei  iiig  an  image  of  a  large  How  lieUI.  this 
te<iiiii(|ne  |)i«mIuc<*s  a  very  comph’x  display  that  maki*s  it  dilfinilt  to  interpret  the 
image.  Finally,  this  te<'hni(|iie  nspiires  the  viewer  to  ’‘lill  in  the  blanks"  to 
interpolate  (low  held  values  hetwis’ii  data  points. 

\\  inkelmann  and  Tsao  were  among  the  first  to  ilevelop  a  raster- ha.se(l 
system  for  flow  fi<’ld  data  display.  Tlndr  approach  took  advantage  of  a  raster- ha.sod 
system's  ability  to  fill  in  and  color  soliti  area.s  on  the  ilisjilay.  This  capability 
alloweil  them  to  show  full  internoflal  value  interpolations  and  proviile  <pialitative 
information  by  associating  sp«’i  ihr  colors  with  (litferent  i>arameter  values. 
W'inkelmann  aiul  Tsao  dc’iiionst rated  the  a«<inacv  of  tln  ir  system  by  comparing 
computer-generated  images  of  (low  (iehl  «lata  with  photographs  from  wiinl  tunnel 
tests;  a  very  closi’  correspondence  between  the  two  image's  iinlicati'd  tin*  valielity  of 
this  techniepK*  (21:177).  Iloweve'r.  tln'ir  system  was  limiterl  to  2l)  displays. 
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pi'orliK'i’rl  a  l<‘('lini«|iif'  tlial  all<*nip(s  l<»  iih<‘  iiiiagrs  siniilar  to  \\  iiikniinaiiii  anil 
i  Mari's  to  liiiild  a  'll)  iiiiag*'.  I'lnMi'  approar  li  unwiali's  “sli«'«*s“  of  flow  liold  imagery. 
tli«*n  ainl  displays  (Immii  in  (lirir  |»»>silioiis  rr’Iativr  to  (lir  How  ln*ld,  as 

shown  in  FiRiin'  J  (  J’JMffT).  AIiIioiikIi  fliis  approai  h  is  iiion'  iisofnl  llian  (In'  “Hat 
iiinif^r's  pi'odni'cd  l»\  W  inkolinann  and  I  sao.  i(  si  ill  Icavi's  >iaps  in  tlin  iinaj>r,  I  Inis 
fort  ini^  (In'  iis«'r  (o  j^in'ss  wha(  (In'  Ihm  Ih'hl  harks  lik<'  widiin  (lirsi'  voiils 


Kioiii  ( 17:12)  I 

Figure  2.  Tlirt'e  Dimnision  Image  Fornn'«|  |)v  Sr'vrral  2D  Images. 

/\  r«)mple(«'l\'  difTf'renl  approai  h  lo  renderin.s,  ilala  was  dexeloped  lie 
William  Fori'iisi'ii  and  llar\e\  ('line  lor  rendeiine,  niedieal  dal  a.  I  Inii  "Maiiliing 
Cilln’s'  algoridim  (Hi)  assumes  die  dala  heiim  used  repiesenls  a  volume  eoiiipiised 
irf  nihie  volumes  ("voxels').  I'^aeh  voxel  has  eit>h(  verdres:  a  speeilie  dala  value  is 
assoeialed  willi  eaeh  verlex.  Once  (he  volume  s  oi  ieiilal  ion  relalive  (o  (he  user  has 
In'en  esiablished.  (he  rendering  process  begins.  K.ach  row  of  voxels  is  examined 
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Ihr  rxistrim' of  lln's«‘  values,  a  is  voiisitlertsi  a  “()'  or  a  “I”,  (l<'|)«‘ii(liiig  on 

wliet  ln*r  ils  <la(a  value  is  below  or  above  ibe  vabie  IxMiig  soiiglil.  Onee  all  t'iglil 
vert  ires  liave  Ims'Ii  |>rores>MHl,  a  siirfar  e  is  <|evel(»|Ms|.  This  is  arroin|>lislietl  by 
|>eiTorniing  a  (able  ltK>kiip  bas«s|  upon  the  patlerii  foriiKsI  by  I  be  vertex 
rla.ssifirat ions.  I  bis  bH>knp  will  redine  tbe  snrfare  lo  oin*  of  fonrlts'ii  possible 
snrfares  that  ran  be  prodnosi  by  linking  ibe  voxel's  edge  midpoints.  After 
assigning  a  clirerti«>nal  normal  base<l  on  tbe  snrfare's  gradient  and  orientation,  a 
lighting  model  is  applied,  tin*  snrfare  is  rendered,  and  pror«*ssing  moves  on  to  tbe 
next  voxel  (KirKil  -  16G).  'I  bis  algorithm  represents  a  railirally  new  tecliniqiie 
kimwn  a.s  volumetrir  remlering.  Its  advantage  is  that  tbe  ilata  iiseil  is  a  3D  set  of 
values  rather  than  a  collertion  of  planar  polygons.  However,  tbe  marrbiiig  cul)es 
algoritbm  operates  by  assuming  that  any  intersecting  planes  will  pass  tbrougb 
either  tbe  edg<*  midpoints  or  tbe  vertices  of  any  given  voxel.  Hasml  on  this 
assumption,  the  algoritbm  then  reverts  to  using  planar  polygons  to  province  a 
sMrfar«>  of  interest . 

One  problem  enronntered  with  31)  terbnirpies  is  that  they  have  traditionally 
tried  to  lit  geometric  primitives  su<b  as  triangles  and  (piadrilaterals  to  tbe  surface 
of  an  object  (3: lUtl- lid).  'I  bis  t<‘clmi<pie  provides  tbe  ailvantagt*  of  well-known 
iiM'tbods  for  deriving  ess«’ntial  surfao*'  information.  However,  there  are  several 
drawbacks  to  this  approa<'b  wIkmi  att(*mpting  to  render  a  flow  fiebl  imagf'. 

I’o  begin  with,  using  this  terbnirpie  can  be  rlilPicult  when  dealing  with 
com|)lex  surfaces;  “boles”  can  appear  in  tbe  smfa<e  being  rendered  due  to 
t ridirect ional  pol\gon  fitting.  Ibis  problem  is  analagtms  lilting  sf|nare  p<'gs  into 
roiiinl  boles;  when  a  polygon  is  use«l  to  build  a  snnM)lb  siirfa<  e  sin  b  as  a  sphere, 
wrinkles  and  gaps  may  appear.  W  ln'rev<M  these  problems  oci  nr,  spe«  iai 
interpolation  ronliin's  need  to  !«•  invoked  in  an  I'ffoil  to  "pali  b"  the  surface 
(l!t:17U).  'Ibis  interpolation  may  b'ad  loan  inaccurate  image'. 
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\  lli»t  lnT  ill  .i\\  li,n  l<  III  |ii  i|\  i;iiii  lit  I  iiii;  !<■(  Iiiiii|ii('v  t  lie  I  ( -si  Ml  I  mil  I  lii'v  |il,ii  r 
oil  till*  \  i<‘\viii);  of  iiii<l('rlyiii|^  iiironiiatioii.  ('iiI(*nk  tiiiir-iiili’iisivi*  sorting  roiitiiK's 
an*  iis<s|,  only  I  In*  oiitrnnosi  slinll  of  data  can  lx*  ilisplaytsi  and  liai  kgroinid  diMails 
am  nii(sk<'«l.  A  liiial  drawliack  to  lln'sr  Its  lnii<|ii(‘s  is  llial  lin'V  tend  to  In*  liascd  «m 
liinary  dtsisions.  I'lnis.  a  portion  of  I  lie  litdtl  t‘itlit'r  is  tir  isn't  rontaiiitsl  williiii  the 
sin  fact';  no  allowainfs  am  niadc  ftir  “'fnz/y  zoiit's"  wlinm  portions  of  lln>  field  am 
appitixinialfly  (Inil  not  exactly)  etpial  tt»  the  isovaliit'  seltsletl  for  rendering.  I  liis 
short  ((lining  nia\  mask  critical  slit'ar  areas  and  dt'privt'  tin*  engiins-r  of  erneial 
iiifurniation.  Herein  lie  the  weaknesses  of  (he  iiiarehiiig  eiilx's  algorilhin:  (1)  the 
desired  data  vahie(s)  will  not  iieeessarily  pa.s8  through  Ihetslge  midpoints,  and  (2) 
enrved  surfaces  may  only  he  npproximaltH  Uy  triangles.  Furthermore,  although 
surfaces  are  common  in  medical  images,  the  only  real  surfaces  within  a  flow  field 
are  those  of  the  object  being  tt'sttsl.  The  flow  field  itself  has  no  surfaces  prr  sr;  It  is 
a  volume  of  scalar  Valin's.  However.  I  Iwlleve  (hat  volumetric  rendt'iing  (echnitpies 
may  be  ii.sed  to  firodnee  accurate  and  usr'fu)  image.s  of  3D  a<'rodyiiainic  flow  field 
data  by  (h'velopiiig  surfaces  of  isometric  values. 

Hobt'ft  Haber  has  stated  that  "'\olumffrir  it ndt  rhitf  methods  are  an 
important  new  development  for  visualizing  three-dinn'iisional  systeins  (  Fidt'i).’' 

Man  bevoy  has  developed  a  volumetric  rendering  techni«|ne  for  medical  aiul 
eh’clron  density  visualization  (I'lritM-  This  technitpie  doesn't  r»'ly  on  gpomt'tric 
primitives;  instead,  the  volmiie  of  data  is  rendered  hy  using  ray-tracing  technitpies 
(as  described  in  (10)  and  (23)).  Simulated  light  rays  are  fired  through  the  array  of 
voxels.  As  a  ray  pa.sses  through  a  voxel,  the  data  value  at  the  point  of  intersection 
lK'tvv«'«'ii  the  ray  and  the  voxel  is  Irirmearlv  interpohiled  from  I  he  dala  valiu's  at  the 
voxel's  eight  vert  id's,  ftased  on  (his  value,  the  (olor  and  1 1  aiislusi  cm  c  of  (he  voxel 
is  de(('rmin<'d  and  fa(  (ored  in(o  (he  rav  s  linal  (olor.  \\  hen  (he  lav  has  passed 
(hioiigh  the  eiidre  volinne.  (he  (olor  of  (he  ray  is  (raiisferred  (o  (he  pixel  (hrough 
which  (he  ray  passes.  When  all  pixels  have  been  (olored.  (he  image  is  romple(e. 
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ln'ing  foi«  r«l  into  a  siufacr.  iii>  ga|»s  will  iw  gdirial«‘<l.  1  lu’  alisniri*  of  gai)s  (or 
|>al«'lii>s  covdiiiR  j(a|>s)  means  that  an  iinaRe's  aernrarv  is  liniitdl  only  l>y  the 
res«>lntion  of  (he  display  device  heinji;  iistsi  and  (he  d(‘nsi(y  of  the  data  mesh. 
Se<oinlly.  (he  ( ranslnsceiice  of  a  vox«'l  is  de(«-rmin('<l  hy  its  value  s  proximity  to  a 
pi(‘sele«  (ei|  target  value;  areas  where  data  values  are  "close”  to  the  target  vain*'  will 
lie  “slighllv  '  colored  (I'itd'i).  This  <'olor  gradation  gives  llic  reseaicher  visual  <  |in's 
ahoiit  t  he  grailient  ainl  t  hi»  kin'ss  of  an  isovalne  surface.  Object  ( ransinscein  e  also 
permits  a  view  of  underlying  and  hai  kgronnd  strncliires.  As  L<*voy  notes,  his 
technirpie  may  he  “easily  iiKnlirn'd  to  render  interstitial  volmnes  as 
semi-transparent  gels  (lOtlh)”.  Ha-sed  on  lievoy's  ix'sults  and  the  fact  that,  the 
eh'ctron  density  data  set  he  ns^xl  is  similar  to  a  flmv  held  data  set.  I  believe  that 
this  (echni(|ne  can  l>e  adapterl  to  How  Held  visualization  to  prodnce  3D  imagr’s  of 
an  <'ntire  a<*rodynamic  flow  field. 

/../  linsir  H< iiiiirr  nirnto 

'lo  be  a  fnlly  effective  engine<’ring  t<K>I,  an  image  generator  should  have  the 
following  capabilit  i(^: 

a.  Data  shouhl  be  displayetl  throughout  (he  entire  flow  field  image  rather  than 
just  in  isolated  areas.  This  capability  provides  a  complete  image,  thus 
assuring  ( he  e?igineer  that  (jotr'iitially  critical  features  are  not  being  omitted. 
Conversi’ly.  the  user  shouhl  be  able  to  isolate  specific  subvolumes  in  order  to 
*‘liminat(‘  regions  of  little  or  no  inteix'st. 

1).  The  objec  t  being  (estc'd  should  be'  inc  ludc'd  in  the  image  along  with  the'  How 
field  data.  The  di.s|)lay  of  critic  al  flow  lic'ld  (eatnres  is  of  little  value  if  (hc're  is 
no  way  to  correlate  these*  feature's  with  (heir  euigin  or  how  (lie'v  alfc'ct  the' 

Ol'T. 
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lx*  |>i'i‘M«*iil  williiii  a  How  If  all  valiios  aio  iliNplavi'tl.  tlx*  iiiiagt*  will  l>e 

c'oiifiisiii^.  ir  l«K)  r<'W  an*  (lis|ilayi*<l.  <  rili«'al  <)ala  may  lx*  omit liy  allowiii;^ 
tlx*  iisri'  to  soIitI  laim*l  valix*s  for  s«*U*<  t«*«l  |>aram(*lf*r.s.  (lx*  systnii  can  r«‘ixln‘ 
imaf^cs  that  isolate  s|x*cilic  valix*s  <»f  int«*f4*st. 

rix*  i'an)><'  of  acci'jit alilt*  ilcv  ialioii  from  a  target  paratnctci'  xalix*  must  In* 
iist'i -<  (lilt  rolli'ii.  I  liis  f(*aliii<*  allows  tlx*  iis<*r  to  focus  tlx*  imaging  syst(*m  as 
fiix'lv  as  lx*  wishes. 


DilTerenl  colors  slionUI  lx*  as.soriat(*(l  with  each  target  value  *»p|eet«l  l*y  the 
user.  Without  color  <lin’«*rentialion,  there  is  no  way  to  distiiiguish  one  target 
value  from  aiiotlx*r.  A  color  k«*y  should  al.so  he  |>rovide<l  in  the  image.  This 
key  should  identify  tlx*  \  alix*  ass<xlate<l  with  each  color  in  tlx*  disiilaj*. 

Helat  ionshii>s  between  |)arameters  should  he  user- con  I  rolled.  I'liere  are  Hve 
ha.s;c  flow  field  parameters:  pressure.  energ\‘.  and  thr»*e  flow  velocity  vector 
valix*s.  All  other  How  (i«*ld  valix*s  may  he  <h*riv(*d  from  these  parameters. 
llow<*ver,  the  user  must  Ik*  able  to  sp^vify  how  these  parameters  are 
int<*rrelaled  lx*fore  tlx*  vahx*s  can  lx*  g«*ix*rated  and  the  results  displayetl. 

riie  syst(*m  must  allow  tlx*  us(*r  to  sp(*cify  tlx*  vi**wing  g(*ometry.  d'his 
involv«*s  positioning  the  viewpoint,  the  viewing  plane,  the  volume  lM*ing 
reixh*r(*d.  and  all  light  sources. 

Temj)oral  features  shoidd  he  displayed.  Many  flow  field  Matures  (such  a,s 
vortices)  (h*v<*lop  and  dissipate  over  a  perxxl  of  time.  In  order  to  study  these 
f<*atiiies,  theii cnliri*  lif<*spaii  must  he  rendered.  Ideallv.  the  ns('i  should  lx- 
ahh*  to  specily  a  start  and  slop  lime,  along  with  ihe  nnmher  of  fram<*s  lo  lx* 
g<*ix'ral(*d  wilhin  lhal  lime  span.  (liven  ihesf*  inputs,  the  sysl(*m  wotdd 
out  pul  I  lx*  r**(piired  nnmher  ol  imag(*s  to  prvxlina*  Ihe  sp(*cified  “movie". 


I.  |{itlii'il  \\<i|||  li.i-  olisi'i \ri|  I  li.il  .  t  •liiii<-n->i<>ii>il  \  nliiiiii-l  I  ii  ...invulvi' 

V  isuali/.iit  ion  lorliiiitiiK'N  llial  do  iiol  (  iirmitiy  (‘xist  in  a  low-cost  workstation 
<Mi\ ironnirnt .  Novcrtlick'ss.  without  nsiiiR  tln*s«'  inrtliocls,  it  is  ilifiicnlt,  if  not 
iiii|u>ssil)l<\  for  a  scientist  to  Kain  a  ^(mmI  nn<l«*rstan<lin|^  of  a  wide  variety  of 
<lata  . . .  In  tnder  to  perniil  nsaJ^e  hy  (lie  largc'sl  possible  nninher  of 

ciij^iin'ers  and  resi'archers.  a  system  should  be  d»‘sij»ne<l  for  use  at  a 
workstation.  It  should  also  be  as  portabh*  a.s  possible  to  allow  its  ns('  on  a 
varii'ty  of  systems. 

/./  .Seopr 

Since  this  thesis  is  an  initial  attempt  to  n.s*'  voinnietric  lechnicpies  to  display 
(low  held  «lata,  I  will  not  be  rieveloping  the  entire  graphics  system  outlined  above. 
Instr’ad.  I  will  dt'sign  and  implement  what  I  consider  the  basic,  most  essential 
featnr<'s  of  the  system. 

rite  first  capability  I  will  provide  is  the  volumetric  display  of  the  How  field 
data  along  with  tin*  objeit  uinh'i  test.  I  will  also  implement  mechanisms  to  allow 
tin-  user  to  specify  the  desire<l  viewing  geonwtry  along  with  any  specific  snbvoinme 
of  interest . 

Another  capability  1  will  providr-  will  allow  the  user  to  select  target  values, 
value  variances,  and  color- value  associations  (se«*  specifications  c.  d.  ami  e  above).  I 
will  also  provide  a  color  key  in  (he  image. 

Due  to  tim<’  constraints,  I  will  not  attempt  to  provide  the  facilities  for 
spe<  i(ying  parameter  interrelationships  or  antoniat ically  g«'m'ra(ing  («'mporal 
leatnre  dis|>lay  se(|uences  (syst«Mn  capabilities  f  and  It.  respec  t ively ).  I  he  first 
(  apability  wonhl  recpiin’  the  tlesign  ami  implementation  of  a  command  language 
t  apabh*  of  allowing  the  nsei  to  both  sc-hs  t  parameters  and  ilest  libe  how  tliev  arc*  to 
be  manipniated  to  derive  a  new  paraim'ter.  Instea«l.  I  will  limit  the  system  to 
displaying  one  scalar  parameter  per  image.  The  second  capability  (time  frame 
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is  llml  I  ilun't  Iw'lirv*'  tlini'  is  siillirirnl  linw  MVAi);<l)|i‘  for  nir  to  ilrvolop  tins 
fr;iliiro.  Sim  oihIIw  oiici'  iin|tlrini'iilr<l.  t  li«'  systi'iii  I  Itiiilil  will  l»»“  alilr  to  gi'iu'iHlo  a 
liinr  s«s(ii<’iir  f*:  liowi'v*'!',  llio  iisi'i  will  lia\«'  l«»  s|MH  ifv  «‘a«  li  frame  iiiiliv  iiliially. 
iVIl  lioiif^li  I  his  is  a  mor<'  tislioiis  iiietliod.  it  will  |iro<liH'(‘  tlie  ilcsireil  results. 

i‘'iiiall\'.  I  will  strive  to  make  tiu'  system  as  portahle  as  possihle.  I'liis  may 
r<s|iiir«'  for«'goiiift  system-speeilii  utilities  in  favor  of  writing  more  geiu'ializeil 
routiiu's;  however,  the  gain  in  support ahility  shoiihl  offset  this  cost. 

/..>  I'hfsis  (Ji-gnnizrttion 

'I’Ik'  r<Mnaiii(ier  of  this  th<‘sis  is  orgaiii/anl  a.s  follows: 

Chapter  2  outlines  the  <le<isi<ins  n‘<jiiir«Hl  to  ilesign  the  system  outlined  in 
S<‘rtion  l.l.  In  this  chapter.  I  discuss  the  alternative  methofis  I  considere<|,  then 
explain  why  I  sele<  te<l  the  options  I  u.se<l. 

('hapter  3  descrihes  the  iiiiph'iiieiitat ion  of  the  imaging  system.  ,\long  with 
a  general  overview  of  the  system  organization.  sp«-cific  key  algorithm  details  are 
discnssi'd. 

Chapter  -!  summarizes  my  tliesis  results  along  with  recommeiulations  for 
future  system  uses  and  enhanci'ments. 

Finally,  a  user's  guide  for  my  system,  the  \’olumetric  Imagery  Program  for 
Fiigineering  Research  (VIPFR)  is  iuchuled  as  an  appendix. 


II.  Si/.striii  Dt  stiji! 


I)i  .lost-pli  Sllall^.  Hii  i«  al  al  llii’  Air  I'cm  c  W  riglit 

Arronaiit ical  LnlMnatorv  ( AKW  AL).  ly  iu)ti‘«l  (hat  “1o  ki'»'|>  up  with  llu-  r«’st 

of  IIh'  \vorl«l  . . .  wu  to  <l«‘v«*l<)p  a  g»MMl  tlir«'<'-<lim*‘iisi()ual  grapliirs  rapaliilily  for 
<lisp|ayiiig  (low  (irIHs  (  Jl)".  A  paiwl  rrport  to  tlw  National  S«  imn’  l‘'uuu«lat ioii 
ri'arlu's  iniK'li  tlir  saino  (oiK'lusiou.  noting  that  “tlinr*  aro  many  known  st  nut  hits 
assoi'iatrd  with  llniil  Mows  sin  li  as  slitx  k  wavi's.  vorlirrs.  shoar  lay<*rs  and  wakes 
. .  .Software  algoril Inns  ami  visnalizalioii  (<*«liiihpi<*s  are  m’orh'tl  to  identify  these 
strnel  nil's  in  ;U)  How  fiehl  solutions  . . .  Innovative  ti'rhnirpies  for  displaying  31) 
data  are  reipiired  l>oth  in  hardware  and  software  (l8;Ar))“.  1  believe  that  a 
relativi'ly  new  volninetrie  remlering  terhnir|ne  pres«*nts  a  solution  to  (his 
visualization  problem. 

Mare  Levoy  ri'eently  ontlimsl  a  teehni<|ne  he  used  (o  generate  images  of 
moleenlar  struct nres  ( 1  "i:2tt-37)  and  noteil  that  most  previous  systems 
appro.\imat«‘d  snifaies  by  arranging  gi'ometri*  primitives  (such  a.s  triangles  or 
rectangh’s)  in  such  a  way  (hat  they  were  “ehrse'*  to  (he  desired  surface.  As  earlier 
noti'd.  this  (<'chni<pie  has  si'vi'ial  major  drawbacks.  N’oinmetric  remh-ring  (VR).  on 
the  other  hand,  is  wi'll  sniti-d  to  (he  computational  fluid  dynamics  (C'I'D)  fiehl. 

/  Dnla  SInirtiin  Sihrtion 

One  data  structure  nseil  in  \  R  is  very  similar  to  (he  data  gem’iatr'd  by  flow 
field  analysis.  I'his  data  structure  is  a  siniph' 31)  array  ("spatial  eiinnierat ion" ) 
consisting  of  data  points,  eai'li  of  which  contains  discn'ie  (low  lii'ld  paraiiK'ter 
valiK's.  I’he  alteriiat  i\e  «lata  strnctnre  for  \  H  use  is  I  he  oilree  (’1:17). 

Oct  rei's  are  t  ree  st  nu  t  nres  t  hat  ha\«'<'iu,hl  chihlreii  per  parent  node,  lo 
biiihl  an  o<  tre<'.  a  rei  nrsi\e  routine  is  used  to  d\ naiiiically  snlnli\  ide  voIiiiik's  into 
octants  until  a  pretlefim'd  level  of  data  homogeneity  within  each  snbvolnine  is 


.11  Ini' \  ('ll  I  I  I  :  I  'i  I.  (  )mi  I-  t  Ills  >.1  Mi<  I  nti-  I-,  limit .  il  m.ix  I  ir  imsiK  iii.iiii|  iiiI;iI  imI  .muI 
tia\*‘rs«>«l  for  iinaKr  ^nirratioii  (!2n:UMI). 

\\  lu'ii  sliort  rxiM  iil  ion  I  inn's  an'  «l«'sin'«l  and  1  .(ttMIJMMI  or  inon'  <la(a  |»oin(s 
ar«'  to  Im'  proct'ssi'd.  rlliricnl  data  inanipulalion  is  a  kny  roiiK'in.  l’|>oii  rradiiiK 
Sandor's  o«lr<'«'  dnsrript  ion  and  al^oritlnns  (20).  I  In'  ot  In'*'  was  my  initial  «|i<ii<(' 
for  tin'  primary  data  slrinlmn  in  iny  sysli'in.  Ilowovor.  I  iijimoto.  Tanaka,  and 
Iwala  liavn  shown  that  I  In- o<  I  ns-'s  clli*  inn  y  j*n»ws  ronsiih'ialdv  worse  than  that  of 
spatial  ('inmn'rat ion  a.s  tin'  depth  of  tin'  volume  heiiif*  lendered  riows  (!(;2I  ).  Siini- 
my  ohjr'ciivo  is  to  render  an  image  of  an  entire  (low  field.  I  antieipatt'il  passing  rays 
Ihrongh  a  volume  with  a  (h'ptli  of  KMt  or  mure.  Oasnl  on  the  fnnlings  of  ITijimoto 
(/  fi/.  I  deridnl  to  ns*'  spatial  emnneration  as  tin*  priin'i|>le  <lata  strin  tnre  for  my 
■system. 

J.  J  Ohjtrt  hijfclion 

Having  s*'h'«  te<l  spatial  ennmeration  for  th*'  syst»'m.  the  n«*.\t  s(«*p  was  to 
s«'h'<  t  tin'  Ix'st  mi'lhod  hit  n'inlering  tln'ohje<  (  mnh'i  test  (Of  1)  in  tin'  satin' 
imag*'  as  the  How  (i*'ld  <lata.  Kurt nnati'ly.  tli*'  ai'rodynamie  ('iigiins'is’  nu'thod  of 
gi'in'rating  fhiw  (i*'ld  «lata  sngg«'sts  a  simph'  nn'thod  hn  i«lrntifying  tin'  Ol' P  within 
tin'  (low  li«'ld.  In  ortler  to  dev<'lop  th«'ir  «lata.  tin'  ('iigiins'is  hnild  a  -(I)  digitization 
of  tin'  Ol'T  ami  use  this  mo<h'l  to  hnild  tin'ir  volnm*'  of  How  (iehl  data  (2).  I  lms. 
hy  using  spatial  ennmeration  as  tin-  hasi*  «lata  strin  tnre  for  my  sysl*'ni.  the 
iin'thud  of  OT'I  repr«*s*'ntat ion  is  also  «'stahlisln'd.  If  tin’  engiin’f’r  who  *l«'\i'lo|)s  tin’ 
•  ligitiz*'*!  moih'l  of  the  ()!' I’  also  uses  a  nni<pn'  data  value  to  ifh'iitifv  tin’  prrs«'n<e 
of  tin'  ()T  r  at  a  <lata  innie.  he  nia\  also  s»'|e«  1  lhal  value  a'>  a  laigel  value  to  he 
sought  ainl  H’lidered  hy  tin’  \  olinin  l  lie  Iniagerv  I’rogram  for  I higineei  ing  Uc'seaieh 
(  \  IIM'.H  ).  I  his  eapahilit  V  will  also  prov  ide  t  he  opt  ion  <  if  rendei  iiig  an  image  i  In' 
flow  field  with  or  withont  tin’  Of  I  heing  di'-piave<l.  It  slninld  he  noli'fl.  howr'vcr. 
that  this  iiK’tlnxI  lepn'seiits  a  d<»nhh’-<’«lg«'«i  svvont.  W  Inie  it  prov  iih's  a  simple  ainl 
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value  UN4>«I  to  |-e|>reseiit  the  Ol'  l'  iiui.st  l>e  iiiii(|ue  within  the  flow  field; 
otherwise,  portions  of  the  How  Held  will  he  rend«‘r*'d  with  the  same  rolor  ns<'d  to 
denote  the  OH'I'  and  an  element  of  amlnKiiity  may  Iw  int rodnr«‘«|  into  the  image. 

J.i  l^nrtiiliHHy 

Since  N'll’KH  is  inteiuh'd  for  eventual  ns«‘  hy  the  «“ngin<*«'ring  community  in 
general,  system  i)ort ability  wa.s  another  «  ritical  <onsi«leration.  To  he  portahh'.  tin* 
system  n<'*'<led  to  he  implemented  with  a  gt'iieral  pni|)ose  high  oialer  language'  and 
as  fr«'«'  as  possible  of  systein-de'pe'inh'nt  f«'atnr»*s.  Once  installed,  t  he  syst<*m  also 
had  to  Im'  as  simple  as  possible  to  maintain.  Since  many  graphics  tools  are  written 
in  the  'i"  programming  language,  1  clerid«'«l  to  implement  m\’  system  in  ‘C’. 
Portability  was  also  eidiance'el  by  using  the  I'tah.rle  graphics  routines.  Tln'se 
eonsist  of  a  pidilic  domain  image  fde  format  ami  a  set  of  tools  to  reposition  and 
combine  images.  I'se  of  t  his  system  provnleel  a  common,  inexpensive,  ami 
pre-cs(  iddislnsl  base  for  spa<  e-elli<  i*’nt  image ‘h'seript  ion. 

In  ordi’r  to  be  portable,  user  interfaces  also  ha»l  to  be  geni'ial  in  nature. 
Although  mou.se-  or  icou-ba.sp<|  interfaces  would  be  simpler  for  the  user,  llu'ir  us«‘ 
wouhl  sT’verely  limit  system  portability.  Instead.  I  chose  siinide  rptest ion/answer 
prompts  for  ust’i  interfaces,  l  ids  approach  rerpdres  only  a  display  device  and  a 
k*'\  board. 

Finally,  all  systc'in-specific  parameter  values  (such  as  the  pixel  dimensions  of 
t  h«‘  displav  dc'vice  being  used)  would  b«‘  isolated  in  a  sc-paratc'  hc'ader  Hie.  Ibis 
Would  provide  a  central  Hie  that  <(juld  be  modified  to  tailor  \  IPl'.H  to  an\ 
particular  system  being  used. 

Having  made  th<’  ma  jor  design  dei  isimis.  I  began  implement  ing  a  system 
that  would  satisly  the  goals  s«'t  forth  m  ('haptei  1.  1  he  next  rhaptei  details  how 
this  program  is  imph'im'utc'd. 


III.  Sjjst.i  in  Implnmutation 


riu'  XoliiimM rir  Imagn  v  l‘l•og^^»ln  for  KiigiiHvriiig  U«‘srarc  li  (\  ll’KU) 
prov'nN's  (lir  rapaliility  for  rnnlniiig  full  ;fl)  iinagos  of  at'Kxlyiiainic  llovv  licM  <lala. 
('Iia|t(n  1  oiilliiK's  thr  prcliiniiiaiy  tirsigii  (l<Tisioiis  dial  weir  r(‘t|iiiic(l  l»«‘lon’ 
pmgiam  impU'inoiital it)ii  Ixt  ii.  I  l':is  <-|iapl<'r  dtMails  how  llw*  a<  (iial  prograin  i.'. 

( oust ni«  t«'<l  and  how  it  p<'ifoniis  its  tasks.  Aft<*r  a  Inirf  o\«m \  irw  of  the  \  i«‘wing 
g«'om«*try  iii\olv('»l  in  volnnift tie  rfiKloriiig.  «*a»ii  of  \  Il’KH's  suhriinctiuns  will  In* 
('xaniiiKHl.  rin*  snhfniutions  are  disnis.s<>d  in  the  same  order  in  whieh  they  oernr 
during  piogram  (‘Neeution. 

■  i.  I  1 7f  whiff  ( h  omi  Ivy 

Before  discussing  the  implenuMitation  d<*(ail.s  of  VIPER,  it  is  necessary  to 
ni\derstand  the  viewing  gwinetry  involv«*»|.  The  viewing  geometry  is  determitierl  hy 
tin*  locations  of  tin-  v  iewer's  eye.  the  viewing  plain*,  the  data  volume,  and  all  light 
sonrc<’s  selecterl  l»y  the  user  (se«'  Figure  3). 

i.l.l  \  it  in  r  's  Fyt .  I  lie  viewer's  eye  is  repres«*nted  hy  a  single  point  in 
space.  Its  location  is  selected  hy  the  user.  .Ml  of  the  rays  that  are  ca«t  to  form  the 
final  image  emanate  from  the  v  iewer's  eye.  The  <lire<  tion  of  the  emanating  rays  is 
determined  hy  the  location  of  the  viewing  plane. 

■i.I.J  y  if  wing  l^lnin .  I  lie  viewing  plane  is  an  array  of  picture  elements 
(pixels).  I  he  size  and  location  of  the  viewing  plane  are  nsi'r  si'lec  t able.  .\  ray  is 
<  ast  fiom  the  eye  through  the  c«’nter  of  eaeh  pi\<-l.  Oin  e  the  ray  is  deline<l.  it  is 
extc’iidefl  as  far  as  necessary  to  (h’termine  whether  or  not  it  intenepts  the  data 
volume. 
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il.il.i  voliiMir  iii;i\  lu'  lliiiiii;lil  it|  ,is  ,i  Ihi\  <il 


siis|u>ii<l«'«l  ill  MpiMT.  riiis  Im>x  is  ('onipi'isi'ii  «)f  xoliiiix' «>I<mii<miIs  (\ox<'Is)  arraiiKisI  in 
a  rrrtaiifi^ulai'.  ;ll)  hm'sU.  I'Im'  KMRlinn  nf  lln*  *lala  voIiiiim'  is  iis«>i  sins  irnHl  by 
('iilt'i'iiiK  (li«‘  <l«'sirr<l  c  xHirdinali's  fur  lli«‘  «'«'iil«'i-|>oin(  of  (Ik*  voIiiiik*.  I''ixiiig  this 
jioiiil  also  (lofiiK's  (Ik*  |ilaiK*s  lliat  form  (In*  fat«*s  of  (Ik*  voImiik*.  In  (liis  itsiksM, 
s|»a(ial  (*nnnK*ra(ion  ads  as  i(s  own  boniKiing  s(riK'(nr('.  a  ^rapliics  (i*<-|ini<|ni*  nsi'd 
(o  iin|)i’o\(*  a^oridim  «*(firi«'my. 


Lij^lit  Source 


Dal  a  VoluiiiP 


\  K’wiiig  I’lane 


rignre  3,  inilial  V  lewing  CK*ome(ry 

■i.I.f  Lifihl  Souros.  d  in*  final  cIcniciU  of  (be  viewing  gcoinr(r\.  Iigb( 
source  loca(  ion.  is  anodier  nsi’i-spev  i(i«'d  lea(nre,  l  or  eacb  ligli(  sonrn*.  (Ik*  nsei 
Mnis(  provi*l<*  (lie  <l<*siref|  local  ion  coordiiiales.  |^a«  li  liglil  source  is  nn»d<*led  as 
poini  dial  «*nn(s  parallel  rays  of  while*  lighl.  I  bis  moeb  l  was  sf*le»  (('d  for  ils 
simplic  ily.  OiIk'i  mo<b*ls  max  prov  ide  a  more*  r«*alis(i(  approxiinal  ion  of  ligbl  s 
behav  ior:  however.  (In-y  also  reepiire*  a  higlK*r  coinpnialioiial  overhead. 


t’. /.  »  liuliiil  S( /«/(.  I  lir  iiii(i<tl  \  K'U  inii  (•niiiH'l  I  \  iii.iv  lie  ;il  fiimi 

t  wo  <liirn**iil  iMTs|M'»  tivfs:  tin*  iisn's  aiul  (Ik*  proniHinmiM's.  Altlioiigli  llM*y  Ri** 
v«M  V  similar,  tlioy  liavr  ilisliiirl  «li(r«*roin  «>.  Tli<’  iisrr's  |>(M'sp«’<-tiv(‘  is  (•ni(pro<l 
al>f>iil  IIh'  virwn  's  ryr:  lunvrvrr.  as  sliowti  in  KiRnn*  .‘1.  I  hi'  piDf^ram  arliially 
rc'volvrs  alxnil  flir  milnr  «»r  tin*  ilala  votniiio. 

lo  lli«‘  iis«-i.  (In'  initial  virwinj'  f»«s»nM'tiy  plan's  (In'  t»rif;in  a(  (In'  virwi'i's 
ryr  ainl  all  manipulations  arr  l»as<'«l  on  (In*  x.  y.  z  axc's  as  shown  in  l''igniT  -1.  I'hc' 
virwiiift  piano  is  coiiti'iod  on  ainl  oitliogonal  (o  (ho  z  axis  at  a  nsrr-specifi<*(l 
(listauco  from  (  he  viewer's  eye.  Likewisr*.  the  data  volume  is  initially  ceuteretl 
orthogonally  on  the  z  axis  at  a  nser-spe<ilied  distance  from  the  eye.  Kiiially,  the 
light  .sonn  es  are  snspende<l  in  sjrace  at  their  n.s<*r-.specified  coor<lina(«*s. 

On  the  other  hand,  the  programmer  se«*s  the  initial  gnnnetry  with  the  data 
voinme's  renter  point  at  (he  X,  W  Z  axes  origins  (see  Figure  Although  the 
viewr'r's  <'ye  is  still  point«'<l  <lirectly  flown  the  Z  axis,  it  is  moved  in  t  he  nrgnlivr  Z 
flirertion  to  establish  the  nser-specifh'fl  flistance  from  the  viewer's  eye  to  the  renter 
of  (he  data  volume.  Once  the  viewer's  eye  position  is  flefined.  the  viewing  plane  is 
Cf'iiteif'd  about  the  Z  axis  at  (he  jnoper  flistance  from  the  eve.  Likewise,  the 
Z-«)rdinate  of  each  liglit  source's  position  is  ailjustefl  rt'lative  to  tin'  eve's  Z-f)rdinate. 

■  i.1.6  I  V<  wiuq  Traiisfoniifilions.  In  imler  tf)  |)rovide  (In'  capability  of 
viewing  all  portions  of  (In'  How  fieltl,  (he  standaifl  graphics  viewing  lransff>rms 
(rotation,  translation,  anti  scaling)  must  In'  available.  Tin'  flifference  between  the 
nsf'i's  ainl  the  programmer's  percept itnis  tif  the  viewing  gisniu'try  arose  from  (he 
imph'im'iitat ion  of  thf'st'  transforms.  I'rom  the  ust'i's  point  of  v  iew,  (he  diihi  rohinn 
is  bf'ing  rotatf'd.  ( ranslat«'»l.  ainl  scah'fl.  Ilowexci.  ( raiislorniiiig  (In'  ilata  vobiiin' 
wtMilfl  ri'<|nirf'  changing  «'a<  h  voxel  s  lo<  a(ion.  .\  (v  pit  al  \')bmn'  may  contain  more 
than  a  million  voxt'Is:  at  IN  mat  In'iiiat  i<  al  operations  (nine  ninlt  i|>licat  ions  and  iiiin' 
afiditions)  per  voxel  (his  approach  wtinhl  n 'finin' a  gn’at  tleal  of  computation  time. 
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I  Im'Ii'Ioii'.  (Ik-  ni(,itiiin  .hkI  (  i  •ni-'l.if  mmi  niMTiif  khi"  .m-  |h-i  ftti  im-'l  i»n  (Ik-  vii-urr 
I’vr  posit  ion.  rHcli  light  source's  Ux-at  ioii.  hihI  racli  pixel  coordinate.  In  order  to 
aclii«'ve  the  r«*sidls  desirerl  hv  the  nwr.  the  ohjts  ls  are  t ransrortiKsI  in  the  opposite 
s<Mise  tt»  the  iiwi'h  «'ntri«*s.  Kor  <*\aniple,  if  the  ns«T  wishi’s  to  rotate  the  ilata 
volume  IV  ahont  the  axis.  \  IPKK  will  rotate  the  viewer's  eye.  all  light  sources, 
and  all  pixels  -1')“  ahont  the  .\  axis.  'I'lins.  the  v<»lnme  remains  fixed  whih'  all  tin- 
«*(hei  comptnients  in  th*’  \  i<*wing  g<'ome(rv  are  moved;  howc‘v«-r.  the  ilittsion  of 
changing  the  volume's  orientation  is  maintained  for  the  nsr'r's  convenience. 

Ihi fortunately,  the  scaling  transforiii  conid  not  he  effecirHl  hy  o|)erating  on 
the  eye  and  pixel  positions;  the  vointne  itself  needed  to  l)c  transformed.  However, 
the  nature  of  the  volume  allowed  me  to  scale  it  without  modifying  each  voxel.  The 
key  to  my  sointion  li«*s  in  the  fact  that,  althongh  spaced  irregularly,  the  voxels' 
locations  are  otherwise  nnifortii;  that  is,  they  form  distinct  planes.  Therefore, 
rather  thatt  storitig  a  distinct  set  of  coor<linat«*s  for  each  voxel.  I  was  able  to  save  all 
.V  ordinates  in  a  i-ditttension  array:  all  )'arid  Z locations  were  stored  likewise.  This 
approac  h  allows  me  to  repieseni  the  locations  of  many  voxels  hy  using  relatively 
few  nnmhers  (r.</..  a  million  voxels  may  l>e  defined  with  as  few  as  ItKJ  nnml)ers). 
Ka«  h  location  or<linate  must  also  l)e  adjnsUxI  to  fit  in  with  the  viewing  geometry: 
this  stej)  involves  snhtracting  half  the  volume's  span  from  each  ordinate  value  (e.g. 
if  the  volume's  A  ordinat«'s  range  from  2  to  (i.  1  wouhl  he  subtracted  from  each 
ordinal*’.  Th*’  resulting  span  wouhl  go  from  -2  to  2  and  Ix’  cent*'red  ahont  =  0.) 
When  the  user  wish«'s  to  change  the  size  of  the  volume,  all  that  is  recpiired  is  a 
simple  multiplication  of  «’ach  ordinate  by  th**  nser-sp*’*  ifie*l  scaling  factor. 

■i.J  Tiinift  c.N 

In  addition  to  spe*  ifyiiig  Ih*’  v  i«’wing  gj’omc’ti  v.  \  ll’l'llt  allows  th*‘  ns«’r  to 
s<'le*'t  larg*'l  \alu<'s  of  inte(<'sl.  the  re*l.  gre<’n.  ainl  blue  (lUIH)  color  \alu*’s  to  he 
a.ssociated  with  *‘ach  larg*’t  value,  th*’  maximum  allowable  flevialion  from  the 

I  ft 


vjilurw.  aiiJ  lln-  niiisiniiiMi  loi  sm |;irc>i  iluil  ri>ni’K|ioiMl  with  llu' 

Ihir<*I.  vahu's.  I'lM'fw*  |>arHin<'t«‘ra  ar«*  M(or«*<l  *w  a  |•«t•or^l:  a  w'paiatr  roronl  is 
iltainlaiiif'il  for  racli  lai'K«’l  valiir  soltM'tctl. 

i.J.I  Tnnjfl  \'rtlii(s.  N'll’I'.H  alUms  srl<‘«tioii  of  up  •«>  live  fargrt  <la(a 
valiirs.  If  111*'  us*T  wislu-s  to  have  llu*  Ol' I  l•^'lul<*n‘ll  in  llu*  iiiiaK«*.  (uu*  of  (lu*  laij»<‘l 
valiu*s  must  lu*  llu*  «lala  vahu*  llial  was  us4'«l  lo  <!i*sipiial«*  llu*  Ol’  1'  in  llu*  voluiiu*. 

■i.J.J  Tavijit  Colors.  A  srparale  lUIU  color  tiiplcl  is  sp**ri(io<l  for  ear'll 
l  arj^rl  value  s«*l«‘rl«*<l  by  I  lie  us«*r.  \\’lu*ii  sp<*<  ifviiig  IKIB  I  ripleis,  I  lie  user  imist  take 
care  to  sele<-t  a  rlifTerent  triplet  for  each  target  value.  Failure  to  do  this  will  result  in 
surfai'es  of  difTereiit  isovalues  being  rendered  with  the  same  color;  image  ambiguity 
will  result.  Furthermore,  the  HCB  triplets  shouhl  lu*  distinct  enough  that  color 
variations  resulting  from  surface  highlighting  will  not  create  a  color  cro.ssover 
(where  surfaces  with  different  initial  colors  appear  to  have  the  same  color).  RGB 
triplets  are  «*xpresse»l  as  sets  of  lhr«*«*  r(*al  numbers  that  may  range  from  0.0  lo  1.0. 

i.J.J  .[lloivnbh  Valuf  Drriation.  I'he  allowable  value  deviation  deliiu’s 
llu*  “window”  of  values  that  will  be  reiulered.  lu  order  lo  provide  visual  clues 
about  the  proximity  of  a  surface's  value  to  the  target  value,  surface  sample's  with 
values  that  are  near  the  center  of  the  value*  wiiulow  will  be  more  opaepu*  than 
samples  that  fall  near  the  wiiulow  «*dges. 

■i.J.-i  Oparifij.  Tlu*  final  user-conI rolleel  parameter  is  the  maximum 
sill fa<<' opacity,  riiis  parameter  may  rang**  from  0.0  to  1.0.  with  liiglu'i  \alues 
de'iiotiiig  a  higher  <legr<'<-  of  <;pacily.  An  opa<  il\  of  0.0  will  |>ri«liu<-  surhucs  that 
are*  enl ir<*ly  lranspar<*nl  and  colorless  whih*  a  vabu*  of  1.0  will  be  eomph'lely 
opaepu*.  .\s  mentioiu'd  pre\  iously.  llu*  final  opae  ily  of  a  surface*  sample*  will  eie*pe>nd 
upeui  its  \alue*'s  pre>.\imil\  lo  a  target  value*.  Thisopaeily  is  e|e'le>rmiue’el  by  the 


I<i||<i\viiit^  ion: 


o/x/f//. =  o/Hir////,, 


^1.0-  — 


oiuf  fdtUfl/*  1 

ulloirahlt  tit  riolioti 


for  <iiiy  tlial  falls  williiii  I lir  siMs  ilusI  \vin(lt>w  of  valiu's.  All  saiM|>l<‘s 

dial  do  not  fall  witliiii  l  lir  valni*  window  will  Im*  assif^nod  an  o|>a('il  v  of  0.0. 


It  should  l)«'  noti'd  that  this  opaiily  fonniila  ililfns  from  that  of  l,«'voy's 
( l''»:d2).  I'hr  ri'ason  for  this  di(f«‘rnic<‘  is  tin*  fart  that  Lrvoy  is  trying  to  rcMidri 
snrfai't's  of  a  givi'ii  thiiknrss  with  thr  following  fortnnla: 


M/wW/.V,, =  I 


I  ^  I  m/m  lorgt  t  MompU  | 

IhirkncAA  |  local  gradinit  \ 


(2) 


Aftrr  ex|M'rimenting  with  IkjIIi  <‘(|natioiis.  I  foninl  that  images  rendereil 
using  Kqnation  2  showe<l  the  same  strurinres  as  images  rendererl  with  C(|nalion  I. 
However,  Fx|nation  2  proilures  a  dimmer  image.  Due  to  the  brighter  images,  along 
with  the  fart  that  it  is  com|>utationally  faster.  VIPF.R  uses  Fc|uatiou  1  to 
determine  surfare  oparities. 


■  I.  I  Ifny  Casting 

As  previously  noterl.  earh  viewing  ray  emanates  from  the  viewer's  eye  and  is 
direrted  through  the  renterpoint  of  a  pi.xel  in  the  viewing  plane.  These  two  points 
determine  the  dire<  tion  of  the  ray.  Onre  the  ray  is  formeil.  it  is  normalized  to  a 
magnitinie  of  unit  length.  I  he  ray  original<*s  at  the  pixel  through  which  it  passes. 

I  h<’  first  step  in  ray  easting  (aft«‘r  the  ray  is  formed)  involves  (h'termining  if 
Ihe  ra\  inlerri’pis  tin'  data  volnnu'.  .\  two  step  process  is  re(|iiired  to  pi’ifonn  this 
rhedv.  I'irst.  the  ray's  point  of  origin  is  t«'st«'*l  to  see  if  it  is  hx  aled  irilhin  the  «lala 
volume.  If  this  is  the  ra.s*’.  Ih<’  ray  is  assumed  to  inleiifpl  lh<'  volume  with  an 
imparl  point  that  roinriih’s  with  Ihe  ray  s  point  of  origin.  If  the  ray  d(M'sn'l 
originate  within  the  volume.  Ihe  n«’.\l  step  is  to  sii'  if  the  ray  ext«*nds  into  the 
volume.  To  perform  this  dierk,  the  distance  from  earh  face  of  the  data  volume  to 
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(In*  ray's  poiiil  of  «irif>iii  is  l»v  a  ilin**  <i«»iial  conipoiM'iil  «>f  (li«'  ray  (for 

iiistaiii'r.  (lii>  «listaim's  from  llir  ray's  urigiii  lo  llir  h-fl  ami  rinlil  fa<'»*s  of  (In*  <lata 
volmnr  arr  ilividotl  l)v  tli«*  ,\-<om|K)m‘ii(  of  ||m‘  ray).  Oner*  all  si.x  fa<r-(o-«*yr 
tlislaiM'f's  are  c  al«  iilal«‘«l.  ••arli  ray  compom'iit  is  mulliplitsl  l*y  om*  o(  tim  tiisiamos 
ami  aiMr«l  (o  flu*  rorrospomliiif^  oyo  positional  or«rmat«*  lo  \  i«*M  a  planar  imparl 
r(M)r(linal('.  If  lliis  point  is  on  onr  of  I  hr  voinmr  s  farrs.  the  «lislam  i'  is  sav«'«l. 

WIk'h  all  six  <lislam«'s  havr  Iksmi  rlK’^kisl.  lh<*  imarrsl  ray-voinmr  impart  point  is 
s«'l('rt«xl  as  a  st  arting  point  for  voluim*  (ravrrsal.  If  no  intrrsr'ction  points  arr 
fonmi.  thr  ray  misses  the  <lata  volniiie  ami  proressing  proreerls  lo  the  next  pixel. 

If  a  voinme  impact  point  is  fonmi,  the  ray  is  incrementally  extendetl 
llirongh  the  data  voinme  with  the  data  being  testerl  at  each  step.  I..evoy's  approach 
to  this  process  involves  “resampling  the  ray  at  A' evenly  spaced  locations  . ,  .along 
the  ray  ( An  early  version  of  V'lI’KR  us<*d  this  approach  with  disap|>ointing 
results  (see  Figure  4).  This  ter-lmupre  can  exhibit  strong  aliasing  artifaiMs  and  may 
omit  entire  voxels.  I  Im'  reason  for  this  lies  in  lln*  fa«  I  that  the  prohahility  of 
linding  all  surfaces  of  interest  tle|K‘nds  on  the  length  chosen  for  the  intervals 
hetw<x’n  steps.  If  a  large  step  is  taken,  many  interstitial  stirfaces  may  he  misserl.  If 
a  smaller  step  is  taken  lo  account  for  this  problem,  the  algorithm  will  experience  a 
linear  increa-s<*  in  execution  time  (e.g..  doubling  the  sampling  fretpiency  will  also 
double  the  time  r<'fpure<l  lo  trav<'rs<'  tin'  \ohum*. ).  I'inally,  it  is  impus.sible  to  a.ssure 
that  all  surfaces  of  interest  will  be  found:  to  flo  so  would  rerptire  taking  an  infinite 
nnmiH’r  of  sam|>l<’s  along  each  ra}‘.  As  an  alternatixe.  I  have  de\’eloped  t  he 
following  algorithm. 

start  th<'  volume  traversal  ))ro«i's>.  the  ia\-\(»lume  impiu  I  point  is  also 
use<|  as  tlu'  tirst  rav-\«)x<'l  iulerse<  I i<»u  piMiit.  1  he  point  where  the  ray  enierg<'s 
from  tin*  \<)X<’I  is  tlu'ii  (alculatefl.  Ome  tlwse  two  points  are  loniid.  tlu'  inl<'i\euing 
data  may  be  l<'sted  against  the  user  specified  target  value's  (as  ilesrribed  leelow  ). 
Ome'  the*  e  iirre'iit  voxe'l  has  be’eii  pioce'sse*el.  the*  ray  me)ve*s  to  the*  m'xt  voxe*l  by 
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Figtire  1.  Aliasing  Pro«lin'c<i  l»y  Sampling  at  Evrn  IntcrvalH. 


assuming  (hat  (In*  nurnit  exit  point  roin<i<les  with  the  entry  point  into  the  next 
voxel  (see  l''ignre  5).  I'liis  process  contiiiuos  until  ♦•ither  the  ray  emerges  from  the 
(lata  volume  or  the  light  available  to  the  ray  is  deph'ted. 

i-  i  Data  Sninpling  anti  Color  (’nh'tiloiion.^ 

Once  a  ray  is  positioned  within  a  voxel,  the  data  lying  along  the  ray's  path 
must  he  tested  against  the  target  values  selected  hy  (he  user,  classified,  and  (if 
appropriat<*)  factored  into  the  ray  s  color.  This  se(|nrnc<'  forms  a  pipeline  as 
des<  rilx'd  hv  Marc  hevoy  in  (I')). 

■  i-i-l  I'lilininr  hilt  rgololioii  of  Ihilo  \  olin>.  The  lirst  step  in  testing  ( In- 
data  is  determining  the  data  \alnes  at  (In-  ra\  -x'oxel  entrv  and  exit  points.  This  is 
ac<omplished  In  following  Levov  s  suggest  ion  (l-'):;Jtt)  and  trilinearly  interpolating 


I'  il^iirr  Hav  Travf'rsnl  i’lirouRli  A«lja<’«*iil  \’ox**ls 

llw  valurs  at  jio'nttJi  haarxl  on  ihr  data  valors  at  the  voxi’l'a  vertices.  Once 

tlo'sr  two  values  have  Ihxmi  d<'leriiiinr<l,  rarli  target  value  can  l)e  l»*ste<l  for 
iiirlnsioii  along  (lie  ray  segment  delhusl  l*v  the  entry  and  exit  (mints. 

There  are  lhr«x‘  (rossihle  rateg«»rM*s  that  a  target  value  may  fall  int«): 

•  It  li<*s  In'tween  the  values  of  tin*  two  end|>oinls. 

•  It  <l«M'sn't  li<’  iM-twjx-n  the  Valin's  of  tin'  two  end|»oints  hut  one  of  the 
end|>oints  lies  within  its  window  of  allowalde  va  lues. 

•  It  diH'sn't  lie  helwix'ii  tin’  \alnes  of  the  two  end|K»inls  and  neither  en<l|K»int 
lies  within  its  wiinlow  of  allowalde  vahu's. 
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II  .1  I  .imi'l  V  .iliH'  III  I  Is  mill  I  111-  l.ii  In  I  ,ili-t;iii  \ .  ii  m<i\  In'  ilisn'Hiinlnl  fm  |  lir 
riirmil  vi»xrl.  Ilowcvrc.  if  it  TmIIs  into  on**  *»f  (In*  *)lli**r  two  CHli'Roim's.  Iln*  iHig**! 
vain**  must  In*  |)ai'aiii«‘l('ri:cf'tl  into  a  |■«•lativ<•  «lislaiii**  tfi-ni  ax  follows: 

If  ill*'  larj^*'(  \aln*‘  falls  lM'lw«'eii  tin*  two  ('iiil|>oiiils. 

^  ^  ^ titrg* I  ^  '  ii/rii  ;if  (  i  ) 

^  1 11^  ij-fl  pi  ^  *  ntr^  pt 

I l<»w*'v*'r.  if  til**  targ**!  valm*  falls  into  tin*  s<*coii<l  *  al**}^*)rv  al>ov«*.  I  will  f**|ual  0  or  1 
(l*'|>*'ii(liiif^  on  whirli  **ii<l|)oiiil  falls  within  lli<*  wiinlow. 

Afi**r  all  larg**!  valii**s  hav**  hf'**n  ('|a.ssiri**(l  and  |>arainp1*‘ii7,e*l.  tliry  are 
sort***l  hv  ascpiiding  valiips  of  /.  Aftpr  thp  sorting  is  romplot**.  «'a*  li  targpt  valiip's 
**fr**ct  upon  tliP  ray's  color  may  hp  ralculat<*d. 

Siiicp  this  t**chni(|UP  iiivolvps  locating  tliP  surface  by  r*'lating  its  valiip  to  the 
values  of  th**  ray  s<'gm«*nt's  <*n<l|)oiiils.  the  algorithm  is  free*!  from  the  constraiiits  of 
sampling  tmints  along  the  segment.  As  a  result.  I  have  eliminat(*d  the  alia-sing 
probinn  I  ciiconnicr*'*!  using  "A  «'v*'nly  spa«'«‘'l  lo«'al i*)ns’*.  Also,  sin*'**  lln*  rav  is 
b«‘ing  moved  along  a  vox«’l  at  a  tim**,  the  maximum  mimlH*r  of  tim**s  the  ray  must 
b*’  advaiic*'*!  and  calculations  nia«le  is  dir«*ctly  r*'lat**d  to  th**  niimb«'r  of  vox*'l8  in 
the  volume  instead  of  th**  *lim«‘nsUms  of  the  vohuu*’  itself.  'I'lius.  hir  a  IDO  x  •'W  x  -lU 
volum*'.  a  ray  will  lussl  a  maximum. of  *>nly  ItM)  extensicuis;  it  mak(*s  no  difTerenre 
whetln'r  the  v*)lnm*’  ha.s  a  span  of  I.tMXt  or  l.OttO.IMK)  units.  On  the  *)ther  hand,  if 
VII’KH  used  evenly  spaced  samples,  increasing  the  volume's  span  by  a  factor  of 
l.Utll)  wouUI  also  increase  the  uumiK'r  of  samples  r**(iuir«*<l  by  a  factor  of  1,000. 

■i.j.J  Sumph  Color  and  Oparihj  Axxniniin  nl .  Determining  I  In-  UtIM 
values  to  lx*  us<'d  for  each  target  \ahu'  is  sim|>l<':  I  In*  us*'r  has  alrc'aclv  ch'lined  llu’ 
H(1H  valiu’s  associated  with  <’ach  larg«’l  valu*'.  l)«‘lermining  lheopa«il\'  of  tin* 
surfac*'  is  little  more  involved.  If  th*’  /-\alu*’  asso<  ial*‘i|  with  llu'  largi’t  value  is 
greater  than  0.  tin*  opacity  of  the  surface  is  e«|ual  to  the  user-specified  opacity 

■ir, 
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I. ill'll  willi  llii'  (.iii:,i'(  \.iltii'  II  /  ti  till-  xint.iii'  i>|>.ii  i1\  Is  liiiMiil  li\ 


Niilisl  il  lit  iiig  riiliii ,  lu.ini  r>*r  ,„<ni>ir  i»  K<i»>ttioii  I.  If  /  =  I  llir  lHrg«*l  valiir  is 
(lisrcgiiiilol;  it  will  lir  ni'ciiiitiliMl  fur  vvlini  lh«'  next  voxel  along  (lie  ray's  |mlli  is 
(Moresseil. 

Siirfan  Iliqhliiihinuj.  'I'lie  next  step  in  tlie  reiidei'iiig  |)i|>eliii(‘ 
involvi's  <lel(‘i mining  llie  sinfaie  liigliliglits  due  It)  the  light  sonrees  sehu  ted  hy  IIh’ 
nsi'i'.  I  use  tin*  IMiong  shading  nioilel  to  rleteiniiin*  surface  highlights.  However, 
hefore  the  I’liong  model  can  he  iis<*d.  tin*  normal  to  *!ie  surface  at  tlu'  point  of  ray 
intersection  must  he  d«'termiue<l.  1  use  the  local  surface  grarlient  to  tletennine  the 
surface  normal.  Accoriling  to  hevoy  (I'inH).  the  local  snrfaie  gradient  can  he 
(h'ti'rmined  by 

j(/(-«'e.Vj+i.  “*•)  -  f{  f‘,.yj-\.:k)) 

7(/('»',..Vr  =*:+i)  -  ./■(•'•i.'/r  ■*-!))  (‘M 

riiis  approach  assigns  the  same  surface  giiulient  to  all  points  within  a  given  voxel. 

\  Il’KH  exti'inls  this  techni(|u<’  aiul  uses  bilinear  inter|)olation  to  a|>proximate  local 
gradiiMits  (ainl  thus  the  surface  normal)  as  follows: 

V/(.s»r/uc(  j.  „ )  ss  ((  /( 

1+  I  •  fj.nur  fnc*  •  -  jiir/arf  )  -  f(r„y,„rj  rtcf  *  'Swr/ocf  )  )' 

(/(  •^sur  flier  •  .Vj+  I  •  ~»urfiicr  )  /(•*  »iir /net  >  .Vj  <- ill  r/ac»  )  )i 

(  ./*(  iiir /'Ilf  •  -Viiir  f'lei  I  *  A  +  I  )  ,/*(  .iti  r/'i.  f  *  //iii  f/'ii'f  '  •*  A  )  )  )  (  ) 


Although  \  ll’KH's  teclini«iue  is  more  comjdicati-d  than  hexov  s.  it  gives  a  mor<’ 
aciurate  value  for  the  siirfan’  normal  at  a  given  point. 

Once  the  surface  normal  is  cahulated.  the  I’hong  shading  mo«le|  may  hi' 
applii'd  to  the  surface  to  (h'termine  surface  highlights. 


•'  l  l  I  iillil  \lli  iniiihiiii.  I.tulif  .il Irmi.iJ KMi  ,t  v|ii'i  i.il 

<ihlioiijo;li  MtiiK*  iiM'llioil  wa.s  lo  imo<I4‘I  iIk*  I'lfcctN  ol  light  passing  tliroiigh 

siirfa«i's  of  (liirnonl  o|)a<  ili«‘s.  an  amiiati'  ami  iralistir  inoilrl  wonhl  !>«'  I'Xtronirly 
roni|>l('\  ami  slow.  .\ls«».  as  Mirliarl  K(*«*l**r  iiofrs.  “il*s  n)m  li  mon*  iinportani  to 
rr|>r«'snit  iiifonnat ion  than  to  g«)  for  nalimn  in  onr  lirld  riu'icfor**.  as  an 

all«'nial ivo  t<»  an  inlii<al<’  iiumIoI  of  light -snrfa<<'  int«Ta<  t ion.  I  «l«“\('lo|)r<|  tin* 
following  ni<‘tlio«l  to  a<<ount  f«)r  <h*<i»'asing  light  Irvris.  Cnat tnniatcci  light  is 
assigm'fl  a  light  h'vrl  valm*  of  1.(1.  .\s  carh  pi.xrl  ray  is  hniiu'd.  it  has  an  availahh* 
light  h  'vri  of  1.0.  Ilowrvor.  as  the  ray  |iass«*s  tliroiigh  surfaces,  the  light  level  is 
<le«reased  with  the  following  fornntia: 

liffhl  It  r<  I  =  litjbt  /( IV  /  X  ( 1.0  —  o/sir//.(/,s„r/,„., ) 

'I'his  technif(ne  is  ha.sed  on  the  premise  that  light  may  only  he  reflected  or 
transmit teil;  ttom'  may  he  ahsorhed  by  a  surface.  Another  simplification  I  made 
assnmes  that  surfaces  only  reflect  from  the  sidt*  that  fac<‘s  the  viewer.  Tints,  if  light 
is  rcllected  from  an  underlying  surface,  all  of  the  r♦'llected  light  will  pass  through 
any  intervening  surfaces.  .\  final  assumption  I  ma<le  is  that  all  surfaces  are 
perfi'ctly  smooth:  no  light  is  siattereil  through  reflection.  Although  th<*se 
a.ssnmptions  grossly  simjilify  the  interaction  between  light  and  physical  objects,  the 
resultant  model  iloi’s  proilme  the  effivi  of  light  atti'iination  without  slow.  comph'X 
computations. 

7  Combining  Samiilt  and  Utig.  I’he  final  ste]>  to  be  performed  involves 
factoring  the  color  of  the  surfai’e  into  the  final  color  ol  the  ray.  lo  do  this.  \  IIM'.H 
simply  multipli<‘s  the  H(JM  valm’s  of  the  suiiacc  by  the  opacity  ol  the  surlace  ami 
l  h<’  available  light.  This  color  is  fhet)  Hthlotl  to  lh<'  ray  s  coloi .  I  hus. 

/fClifi  ni  =  IfCliii,,,!  +  l.iglil  l.(  ri  I  X  lf(  I'lis,,,  f  X  (}iniril g^,,,  j„,  , 
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iHmI  till'  lllltll  Ittitil  mI  ill!'  |>|\i'|  lIlMlIll^ll  wimll  till*  |(|\  |l«|>si*si  is  UlM'Il  li\ 
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hrti,  x  x  f  ////>, ). 

•  .=n 

liv  «l<'<  r»'asinj>  ihr  li^lit  |<'v*'l  niul  iii<'t»r|u>rHl  in^  *'Hrli  suiTarr's  color 
coni riUiil ion  to  llir  |ii\rl  as  tin'  ray  is  lrac«'«l  lliroiif^li  llic  volninr.  I  am  able  to 
Ix'iTorm  a  IronMo-ltack  simnnation  ainl  llins  avoi«l  tin*  rxlra  lime  and  o\.'rln'a«l 
ri'<|nir<'d  l)V  I  lie  fronrio-liack  Iravi'isal  followed  l>y  a  liack-lo-fronl  siimmal  ion 
le<  lini<|ne  d«'\<'lo|M'd  l>y  I  riedi'i.  (lordon.  ainl  Ileynolds  (M). 


■/.  )  Siniimni  !/  of  Alfionthm 

•  Head  in  Data  and  Build  the  Data  Volnine. 

•  lnl<‘rarl ively  reail  in  I'ser-Conirolhxt  I’arann'lers. 

•  Transform  the  \'iewing  (leonietry. 

•  For  Eacli  Pixel: 

-  Form  Hay. 

-  ('Iie«k  for  Hay-\olnme  lnl«*rs<*<  tions.  If  Hay  Intersects  the  \oInme, 
I'raii’  tin'  Hay  I'lirongli  tin*  Volume.  At  each  Voxel; 

♦  Lo<  ate  Hay's  F-ntrv  and  F.xil  Points.  Trilinearly  Intevtmlate  their 
V  aim's. 

♦  Locate  all  Hay-Target  V  alue  Surface  Intersections. 

♦  Calculate  Each  Surfaci'  s  Normal.  Shading  Fai  lor.  Coni rihnt ion  to 
Pixel  s  Final  H(1B  Value,  ainl  Light  Vl  tennal  ion. 

♦  Extend  Hay  to  Next  V<»x«'L 

-  Sa\r'  Pixel  s  Eiiial  H(JB  Values  in  a  I  rame  Mulfer. 

•  Sav«'  Frame  Hutfer  in  I'lah.rle  Format. 


IV.  CoucUisums 


I  hr  liisl  llini'  i  liHptrrs  *>f  this  llirsis  |>n*sriil  I  hr  ralioiiair  hikI 
iinpirinriitatioii  <lrlails  inv()l\«‘(l  in  N'll'KK's  ilrvrlopinrnt.  In  this  rhaptrr.  I  will 
assrss  \  Il'MU's  rapahilit  irs  and  prrhn  tnanrr  rhaia<  lr!  isl  irs,  |  will  also  onl  liix' 
what  I  l»rlir\r  would  l>r  maj<»i'  rnhanr<*inrnls  fo  \  IIM'dI  as  well  as  ollirr  Tudds  in 
whif  h  \  IIM'IU  inij>h(  provr  ns«‘fnl. 

{.I  f*fiforinnnn  ( 'hnrarit 

i-l.l  H-rfM  rl  Opinion.  Ina.siiinrh  as  I  am  not  an  arro«lynamirisl  and 
lhri«'forr  not  (inalifird  to  assess  \  IPEH's  utility  in  that  firld,  I  sought  the  opinions 
of  those  who  would  know  what  is  us«-ful  and  xvhat  is  not.  Arconlingly.  I  presented 
N’ll’KU-geiu'iatrd  images  <»f  th(X)rrtiral  (low  field  data  It)  (’apt  Heiau  (an  instructor 
in  tiu'  AFl  r  l)«'partm<‘nt  of  A<Monaut  ical  Kiigiiux-ring)  aiul  Dr.  Webster,  an 
arroclvnamicist  from  tin*  Air  Force  Wright  .Aeronautical  Lahoratory.  Both  of  these 
review»*rs  made  favorable  comments  on  \  IPFH's  imag(*ry  and  capabilities. 

In  an  elfort  to  demonstrate  \  IPKH's  applicability  to  other  fields,  I  used 
\  IIM'.H  to  prvxince  a  series  of  images  ba.s«’<|  on  the  atomic  energy  levels 
surrounding  compresse<l  cry  st  alline  st  met  ures.  Lt  Col  bnpo.  a  professor  from 
AFi  r's  Physics  Department,  was  favorably  impressevi  by  these  images,  noting  that 
they  closely  corresponded  with  the  Jl)  contour  plots  he  had  ih'velopevl  for  the  data, 
lie  flirt  heriuore  stated  that  the  dl)  images  d«’v«doped  l>y  \  lI’ER  were  “much 
belti'r  than  the  polygonal- based  imagi’s  he  had  been  using. 

Finally.  \  IPF.R  was  nsisl  by  ('apt  (laiy  boriinoi  (a  fellow  .Ah  11  student  )  in 
the  ionise  o(  his  thesis  elforts.  \\  hen  he  noted  errors  generated  by  a  Poliiemns 
Navigation  S<  ien<  es  .'l-SP.ACK  ti  acker,  he  took  a  seiies  of  readings  at  speciiic 
spatial  intervals,  t’sing  these  readings.  \  llM.U  was  able  to  generate  an  image  of 
piedicteil  error  values  within  the'  space  where  the  Iraiker  was  be'ing  us«'d.  Studying 


I  ll  <■  iiii.itii's  iilliiwi'il  (  l.iiiiiiiui  1(1  uIciiIiIk  •'itci  iIk  ic^toiis  ()|  liolli  .ind  low 

(•nor  zones. 

Unsed  on  these  a.ssessmenis.  I  believe  that  the  a|)|)lieal ions  of  volnmetiit 
rendeiint!,  slioidd  l)e  fnitliei  explored,  not  only  in  the  held  of  ronipntat ional  llnid 
d>iiamies  hnt  in  other  lields  of  study  as  well. 

/./.J  Aliiorilhnt  nrii.  \  Il’KK  is  a  graphics  program  that  uses  a 

ra\  -(  asliin'  alj!,orithin  to  prov  ide  fidl.  dl)  iina|>es  of  a  vohnneof  data.  In  order  to 
assess  algorithm  elliciency.  one  must  normally  specify  whether  efiiriency  is  being 
expr(^s.sed  in  terms  of  time  or  space  complexity.  In  a,ssessing  VIPPdI.  I  have  found 
that  this  distinction  is  not  necessary:  the  algorithm's  time  nnrf  space  complexities 
are  both  dependent  upon  the  size  of  the  data  volume  (in  terms  of  voxel 
dimensions)  to  be  rendered. 

If  one  disregards  the  process  (»f  reading  in  the  data  volume.  VIPP'H  operates 
with  an  elliciency  of  Ofitr).  where  /)  is  the  nnmber  of  pixels  throngh  which  rays  are 
being  cast  and  r  is  the  volume's  largest  dimension  (in  terms  of  voxels  per  edge). 

I  In*  elliciency  notation  f(»r  the  data-reading  portion  of  \  IPPH  depends  on  the  unit 
(»f  nH'asnrement  to  be  u.sed.  If  measnretl  in  terms  «»f  total  number  of  mesh  nodes  in 
the  data  vvdnme.  this  segment  operates  at  a  linear  ellic  ien<  y  level  of  Ofd)  (where  d 
is  the  mimbei  of  data  points  to  be  read  in).  On  the  other  hand,  if  one  ns«'s  tin* 
volume  s  dimeiisicnis  (measured  in  voxels  per  r'vige)  as  the  standard  of  measure. 

\  IPP-irs  data-reading  performance  may  bo  expressed  as  (){ijl,  )  where  are 

t  he  .  Z  dimensions  of  I  he  dat  a  volume 

j.I.J.I  l:'.miilinn  linn.  It  is  dilli( nil  to  give  a  meaiiingfnl  ligiire 
for  \  IPKH's  exe(  lit  ion  t  ime  due  to  t  he  number  of  fa(  tors  ( < number  of  v(»\els  in 
(he  volume,  pixel  dimensions  of  final  image,  objefl  magnili(  at  ion.  number  of 
surfaces  being  sought,  surfarc  opa(  ily.  <•(<  .)  involved.  However,  to  prodmc  a 
demonstration  video  that  is  lliirty  sevonds  in  length  re(|uired  (he  rendering  of  !)(K) 
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iin.iU'  "'  I  It'-  \<i|i|iui'  I I’lM l< 'I <•! I  li,iil  viiM'l  III  IfMI  <  (is  X  (ix  Ini'  ii  li)l;il  iil 

IliLMdO  iiirsli  iio<I«‘m;  tlir  \  itUni  sliowttl  I  hr  vuliiiiir  hI  ilifFrinit  Hiigirs  mm<I 

•  liirririil  (iparily  srlliiigs  for  tlir  larg«'»  siirfanis.  Karli  iinag<‘  inrasnr<Ml  x  •'»12 
tlir  Digital  K<|iii|iiiiriil  Corporation  (<P\  Mirrovax  III  grapliirs 
systrins  lorati'il  at  AKI’I  ,  tlir  tiinr  rr«iiiir«Ml  prr  iiiiagr  raiigrd  from  27  iiiiiinti's  to 
approximatrlv  hours.  .\lt hough  thrsr  ligiirrs  obviously  pri'cliulr  tlir  iisr  of 
\  IIM'.H  for  iiitrrartivr  imagr  griKMation  ainl  virwiiig.  tliry  arr  iu>l  prohihilivriy 
largr;  21)  iiiiagrs  ii.srd  to  druioiist rat<*  graphics  capahilitirs  at  N.XSA's  haugiry 
|{*‘srHrch  Cruti'i  rerpurr  ''only  t  wo  hours  of  criitral  processing  tiiur  ( I7:r»)“  on  a 
\  PS-.12  snprrcoinpntrr.  Therefore.  I  Iwlievethat  VIl’EH  ileiuonst rales  a  useful 
image  |>ro<luction  rate. 

J  M(  inory  l{(quhrin(  nift.  Memory  rerpiirement  s,  like 
exei  iition  times,  are  difficult  to  (piantify;  therefore.  I  shall  again  resort  to  the 
iiuMiiory  rerpiiremeiits  I  ohsr'ived  during  the  production  of  VFPER  s  demonstration 
viihx).  riie  storage  r<'<pureuient  for  tiu'  «lata  file  was  lt).7ll  megain  tes  (Mh)  of 
oiiliiK*  iiKMiiory;  however,  it  should  Ik*  not«’d  that  the  data  file  was  storerl  in  ASCII 
format.  While  less  memory  would  he  re<|uirrd  if  data  was  storeil  in  a  binary  format 
this  would  indini-  an  elriiiriit  of  system  dependency  and  thus  limit  the  portability 
of  the  program.  W  hile  execiit ing.  \  IPER  rerpiired  2.117Mb;  this  recpiirement  will 
change  in  proportion  to  the  size  of  the  data  volume. 

(.1.1  finia  Di.^plny.  As  noted  in  Chapter  I.  most  ar'rodynamic  images  to 
«late  have  Ikhmi  restrii  ted  to  eitln’r  2D  ilisplays,  2D  images  su|)erimposed  upon  each 
other  in  a  ID  arraiigeiiK'iit .  <<dor-<  ofle«l  r»’gioiis  on  the  surface  of  the 
objc'ct -under-test  (Ol'  l  ).  or  a  ID  clispla\  of  \<m  tors  or  partic  h'  jjallis.  f  igure  (i 
shows  a  typical  21)  rendering  of  the-  How  licHd  around  a  friclioiiless  sphc’re.  in 
contrast,  f’ignrc'  7  shows  a  stereo  image'  of  the  '■ame  How  liehl  as  rendc'ic'd  l>.v 
\  IPER.  file  different  slic'lls  displayed  in  the*  image  corresjioiid  to  various  values  of 
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I  lie  V  <  ll  M  ll  I  I  II  >1  I  III  1,1  I  ,1  I  M  lilt  I  III'  Nlllil  'I  I  I  III  'Ill'll"  ill  I'  I  Mil  \  III  III  I  'll  ill  I  IV  I  III '  lull  IS, 


'<1  nIiovvii  ill  <».  llow«‘\i'i.  l  ari’fiil  «'oiii|iAriwiii  ImMvvii*ii  IIhim*  (wii 

liKiiri'N  sluiwM  iIinI  \  IIM'lll'fi  iiiiHRt*  in  vvlini  is*  jiiiMlkUMl  in  ( lil).  AIm>  m>lr 

llu*  I'oUiml  iiiiiiiiN'i'x  ill  I  In'  ln\vi*r  IrfI  l  oi  iirr  of  I  hr  iiiiAj^r;  lliry  form  I  In*  iniAgr 
l  ohir  ki'v  iiinl  iinlii  Mli'  I  In*  laigi'1  VHln«*s  Immiir  r«*inlrr«'<l  as  iv«*ll  as  tin*  rolor  with 
vvhirli  llirv  an*  IwiiiR  ilisplayril.  Tin*  pah*  nn*v  s|iln*i'<*  in  lln*  iriilrr  of  (In*  ilisplay 
ri*|ni*si>nts  (In*  ()!' T, 


I'ignri*  n.  Trsllwsik  Imaj’i*  of  f  low  Aoonl  a  Spln'ir 

\  IIM'iU  s  aiiiiiacv  may  also  In*  as'irssnl  liy  <  oiii|iariiip,  (In'  roii(oiir  |)lo(  ol  a 
\vak<*  vorirx  (I'ignri*  f<)  willi  (In*  \  ll’I^ITJ'^•ln•^a(l'1l  iinagi*  of  (In*  sann*  <la(a  sr( 

(l  igiin*  (I.  Tin*  raiifti*  of  \  IPF.R’s  (’a|>aliili(irs  is  simj^rsti'il  liy  I'igiin*  10  wliii  h 
shows  (he*  sann*  virw  of  (hr*  vortrx  r«*inh’iH«l  wi(h  ililTm'Ut  opacity  valin*s. 


P'igiiro  7.  (tnagr  of  Flow  Alwiit  a  Spiiore 

1  lirlievf*  ( lia(  I  «'xaini»l<*s  ami  Um*  pirvioiisly-cilrd  road  ion  from 
l>ot(‘iilial  users  deinonst  rate  \  Il’FH's  |ioif'nlial  for  improving  tlie  stale  of 
a<‘rodynamir  How  field  imagery. 

/./.  /  Aliasiinj  ('.ffu  ls.  As  I  noted  in  Cliapler  3.  tlie  original  st.rateg\'  of 
sampling  at  even  steps  along  ea«  li  ray  proilmed  very  pronounretl  aliasing  artifacts 
due  to  tlie  resultant,  gaps  in  coverage  that  occur  betwecMi  sampling  points.  Locating 
surfaces  a,s  a  function  of  s**gment  endpoint  values  eliminated  this  problem. 

Ilowexer.  there  ai'e  still  sources  of  aliasing  that  shouhl  Im-  a<hlressed. 

lo  begin  with.  IIh’  number  of  voxels  in  the  data  volume  itself  will  have  an 
I'lfect  on  the  degree  of  spatial  aliasing  apparent  in  the  linal  image.  Specilicallv.  a 
given  volume  will  exhibit  f<-wer  aliasing  arlifa<  Is  if  it  is  reprc'sentc'd  by  IDO.tMM) 
voxels  than  it  would  if  it  is  comprised  <»f  onlv  ItKI  voxc'ls:  this  is  another  form  of 


Kigur«*  8.  Coiitour  I’lol  of  \\ak«*  V'ortox. 


nliitsing  as  a  fiiiK  tion  of  sain|>liiig  fn'<|unny.  Ilowfvrr.  I  he  niiiiil>er  of  voxels  in  the 
volume  will  l>e  dictated  liy  the  user’s  ihxxIs  and  system  capahilities:  this  is  not  a 
proMem  that  «  aii  Im*  solved  within  \  IPKR. 

Another  aliasing  eir*^^  that  Naries  with  the  size  of  the  voxels  Inung  renderetl 
is  also  inllneiiccsl  hy  tin*  st<s’pin*ss  of  the  hical  grailieni  surrounding  a  given  voxel, 
the  size  of  the  value  window  Ixuiig  us«*d.  ami  the  opacity  of  the  target  surface, 
(liven  the  right  circumstances,  the  entire  voxel  will  In*  coloreil  and  displayrsi  as  ,i 
rectangular  hlock  within  the  fmal  image.  This  effect  is  most  likely  to  occur  when 
using  a  small  valu<’  window  and  high  opacity  c«K-jliriciil  to  render  a  region 
<»)mpris«sl  of  relatively  large  vox<ds  containing  a  slis'p  local  gradient.  The  usi'i  ina\ 
<lampen  this  form  of  aliasing  liv  spe<  ifving  fairly  wide  value  windows  and  low 
opacity  coeflicii'iits.  particularly  when  s«*tling  these  values  for  n'lnh'iiiig  the  OC'I'. 
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Figure  9.  VlPER-CJcucvatotl  lmag<*  of  WVke  Voviex 


AnotlKT  sourer  of  aliH.siiig  is  |>ro<hi(-i'<l  by  i)oiiil  sHiupling  and  is  dr|><’iKlriil 
U|)oii  t  he  resoliiliou  of  tlir  viewing  plane  (inea-snrrd  in  pixel  diniensions).  Tlie 
•  oarser  the  vi<'wing  plane,  tin*  more  prononiKXHl  the  aliasing  (wliieli  manift'sis  itself 
in  tin*  form  of  jagg«*d  edges)  will  be.  I’liis  is  due  to  the  fact  llial  the  rolor  of  an 
entire  area  (the  pix(’l)  is  determined  by  the  eolor  of  an  inlinitely  small  area  (the 
cross  section  of  the  pixel  ray).  Many  schejues  have  been  advance«l  in  an  effort  to 
remedy  this  common  graphics  probh'm.  Typical  sohitions  include  snpersainpiing 
(firing  more  than  one  ray  per  pixel  and  averaging  the  resultant  IKIM  triplets), 
filtering,  and  "binning"  tli«'  imag**  l»v  averaging  the  <H>)ors  of  ailjacent  pixels.  All  of 
tln’se  approaches  im  nr  a  higher  ( (unpntat ioiial  overhead  and  t<Mid  to  degiade  tin- 
sharpness  of  the  (inal  image.  1  exp«'rim<Mited  with  snpeisampling  ami  found  that, 
although  a  lim'ar  imn-ase  in  prtxcssing  time  was  rec(niietl.  little  impiovement  was 


I'^igurr  10,  U'ako  Vorl<*x  Witli  DifTcrcut  Opacity  Values. 


A  final  sonic*’  of  aliasing  is  identiiuHl  l>v  Drebin.  C'ai|)<’nU’r.  anti  llanralian 
wlitj  note  that  “when  rays  diverge  they  may  not  sample  atijaceni  pixels  (5:72)" 
(Since  a  ray  is  fired  through  trfiy  pixel,  I  believe  this  is  an  error  and  /n>rf  shonitl 
re  ad  rnril.).  However,  giving  (he  user  the  ability  to  adjust  the  the  locations  of  the 
viewing  plane  anti  tlata  volume  provitles  (he  means  to  minimize  this  problem.  As 
the  distance  between  the  viewer's  eve  anti  the  viewing  plane  increa.ses.  the  angle  of 
tlivergence  between  rays  tlecreases.  As  the  tlistance  approaches  inrmity,  the  ray 
tlivergence  approaches  zert).  As  an  example  of  this  technif|ne,  when  making  my 
tleiiionst rat itni  vithst,  1  used  a  'il'i  x  512  |»ix«’l  image  plain- s<’t  at  a  <listance  of  5tMM( 
units  from  the  viewer  s  eye.  (liven  this  < onlignrat ion.  the  maxiiuniu  divergent e 
between  neighbt)ring  rays  is  only  O.OItiJ'’.  If  less  divergence  is  desired,  the  image 
plane  an*l  data  volume  mav  be  moved  to  an  even  greater  tlistance. 


Kigiirf*  II.  Suhvoliiine  IsolHlioti  Via  (’ulling  Plane  Method. 

/./..>  I'aroiiii  li  r  (  onlrnl.  I  lien*  are  a  iiiiiiiher  of  paraineteiH  that  are 
ly  <'on(  rolh'd  l»y  IIk'  user.  Tlir'se  fall  iiilo  one  of  two  g«*neral  ralegorH*s: 
scene  specificat ion  an«l  target  valne/j>ro|MTfy  definition.  After  reading  in  llie  data 
volnmi’.  \  IPKIf  prompts  tin*  user  to  enter  values  for  each  control  parameter. 

f.l.l.l  SrdH  SjH  cijicalion.  This  first  category  of  param«*tt*rs  allows 
the  user  to  ‘‘st't  the  scene".  The  u.ser  must  first  specify  the  viewing  geometry.  The 
parameters  that  control  the  initial  g«s)metry  are  the  distance  from  the  viewer's  eye 
It)  the  (('liter  of  the  data  volume,  the  dislainc  from  the  viewer  s  eye  to  tin'  viewing 
plane,  and  the  location  of  all  desired  light  soiirtcs. 

d  h('  user  may  also  selec  tively  "cut  avvav  '  portions  of  the  data  voluiiK'  in  two 
ways,  d'his  permits  the  user  to  isolate  a  specilie  siihv olume  of  interest.  The  liist 
method  of  volume  elimiiiat ion  is  aecomplislu'd  hy  placing  the  viewing  plane  within 

;i7 


I  lie  \  <  iliiiiH'.  .!•'  'll!  i\\  II  III  I'  iu,iin-  M  I U  I  Mini;  (ir-  <-l  i  .iv  liom  i  In-  v  icw  iiih  |  iI.hh' 
iii.sl«'H<l  III  IliM  \  il'WiM  s  rvi’.  I  was  alilf  to  liaiisrunii  (in*  \  irwiiiK  l>liiii«‘  a  ciMliiiK 
plaiir.  I'lir  sitoikI  inriinxi  of  rliniinal  iiiK  |>oi  tioiis  <if  I  hr  volnnic  involv  es 
selei  livi-lv  mnoviiif;  iiser-spei  ilii'd  slalis  fioiii  any  for  all)  of  (lie  volume's  six  facrs. 
I'iK'iic*  12  shows  the  same  Ihiw  held  as  I'igiiie  7:  however.  llVX  of  the  flow  field  ha.s 
Imxmi  seleel  ivi'ly  removed  hv  lakinj;  slabs  off  of  I  he  froiil  and  hark  of  the  data 
volume.  I  he  renlral  'v'/  of  the  volume  is  thus  rev<‘ale<|  for  examination. 

Once  the  suhviiliimi'  has  Ixxmi  delined.  the  user  must  in|)nt  sealing  failors. 
angles  of  rotation,  aiul  translational  distanees.  I'itially.  the  user  is  given  the 
opportunity  to  rhang<'  the  si/e  of  the  viewing  plane.  Onre  ent«*red.  thi'se 
paraiueti'is  deline  the  liual  viewing  geometry.  AnotluT  variable  to  be  entererl 
speeifh's  the  ambient  light  level  within  the  scene.  Finally,  an  IKIlf  color  triplet 
must  bc'  specified  for  the  image  matt<\  I’liis  is  a  color  to  be  applic'd  to  any  pixel  ray 
that  missc’s  the  data  volume.  Kays  that  do  intcTsect  the  volume  are  initially  black 
(K({H  =  (t.O.tt)  and  change  as  they  enccnuiter  targc't  surface's.  I'hc'  objective  of  this 
feature  is  to  display  ihc'  silhouc'lle  of  the  data  volume  against  a  colored  mat;  this 
allows  the  user  to  orient  the  rendered  surfaces  vvithiu  the  volume  itself.  If  the  user 
wishc’s  to  disable  this  feature*,  he  only  ne«*ds  to  select  a  mat  Klfli  triplet  of  (U.O.U). 

/./.  ». 2  I'diijtl  \aluf.s.  I  hc’ ciirrc'iit  version  of  VI PKK  allows  t he 
usc'f  to  s|jecily  up  to  live  target  value's  to  be  sought  out  and  re'iidc'red.  For  each 
targe't  value*,  the*  user  must  spe*<  ifv  the  value*  itse'lf.  a  winelow  e)f  acceptable  value's 
which  is  (e*ntere*d  about  the*  target  value*,  the*  K(JU  cedor  triple*t  to  be  a.ssoe'iate*el 
with  the*  target  value*,  and  the*  opacity  of  the*  laigel  value*  s  suilace*.  1  hc'si*  value's 
are*  assiguc'd  iue|epe*ne|eut Iv  for  e*ach  laigc'l  value*  to  be*  sought.  llovveve*r.  the* 
re’fle'ct  iv  ily  values  (diffuse*  and  spe*cular  re'lle'e  t  ion  eoe*llii  ieui s )  are*  spe'cilie*d  onc  e* 
aiiel  applie*d  uniformly  to  all  re'iulered  surface's. 


I .  I  it  I'lti  hiltilihi.  SiiM  r  \  1 1  ’  I'.U  I''  Idt  s  w  till  \\  ill  li.ni'  ;i  w  i<l<' 

VHi'irly  of  «'oin|)iil iii^  fAcililirs  hI  iIk'IT  iliH|Mtsal.  i  a  ^n-al  <lral  of  (MiipliHais 

•HI  niaiiitaiiiing  |>orlal>ililv.  A<((H<liiiRly.  I  ••srlu*\vr«l  sys(«>ni-H|M*rifi«  in  fav«>r 

of  using  gnn'rally  availalilr  rnin  lions  (siirli  as  sim*  ainl  ('osiin*  fiini  tions)  and 
siraiglil forward  ‘C  (  tMl*'.  As  |>r^‘viunsly  iiuMitiontsI,  I  In'  final  image  is  slor«'<l  in  (In' 
l'l•dl.rl<■  r<>rmat.  Tin'  use  of  lids  rormal  was  jnsliln'tl  lieeansi'  il  is  a  |>nl>lie  •lomain 
roriiial  dial  is  widi'ly  sn|>|nHl<'il  f«)r  many  differi'nl  sysh'ins.  Allliongli  d«'V«'io|><’ii  on 
a  DI'T'  (JPX  Mierovax  ill  sysli'in.  I  have  Ins'ii  ald<‘  (o  (ransfer.  n*«-om|>ile.  and  run 
\  Il’KU  on  a  Siliron  (Iraphies  :J|.JU  Iris  system  as  well  as  a  \  AX  1 1/785  rom|)nl«'r. 

T-  SiiggtshH  I'.nhanmuf  uh 

\  IPKR  is  nol  inli'inh'd  as  an  nitimale  graphics  (iM>l  fur  engiin't'is.  It  is  also 
in»t  inlendefi  to  he  used  for  (|nan(ilative  measnremenis.  However.  1  do  believe  it 
has  die  potential  to  In'eome  a  vi'ry  powerful  ongiin*ering  tool  by  providing  3D 
images  of  aerodynamir  flow  fiehl  dal  a.  1'hese  images  will  make  rpialitative 
assessiiK'iil  of  How  li«'lds  r<'la(i\<'ly  simph'  and  inlni(i\<'  for  engiin'i'is.  in  (In*  eoursi' 
of  (his  thesis.  1  have  (h'veloped  a  prolotype  for  V1I*1’2H  that  is  simple  to  use  and 
provid«'s  many  <if  the  f«'alm«‘s  an  engiin'er  is  likely  to  r«'(|iiire  or  desire.  Although 
useful  ill  its  eurreni  form,  flier*'  are  sev«'ral  features  that  conhl  be  ineorporatetl  into 
\  Il’KH  lo  make  it  «'v*'n  mor<'  robust  ami  powerful. 

f.J.I  Cirid  Configuration.  VIPKH  operates  on  (he  assumption  that  the 
•  lata  volume  will  be  represent'd  by  a  reetaugular  griil.  It  was  not  until  the  la(«*r 
stages  of  my  (In'sis  work  (hat  I  dis«-ov<'n'<l  dial  a<'rod\ nami*  isis  rari'ly  us«'  a 
r('<'( angular  gri<l:  iiistea<l.  (In'.v  rouliin'ly  us«'  a  grid  (hal  is  in'arly  orlhogoiial  bill 
not  r<'(langular  (s<'e  |■'igur«’  13).  lo  •omp*'Usat'  for  lids  dilh'ieiu’e  in  gri<ls.  1  would 
sugg<'sl  devc'loping  ail  inl«'run'diar\  program  dial  •  •>uld  map  a  !•'•  I  angular  grid  into 
a  giv<'ii  inHi-r<'<  laiigular  giul.  1 1  iliin'arly  iiit'i  polat'  the  ilala  \ahn'  at  dial  point, 
and  produce  a  new  data  file  that  routaiiis  a  reetaugular  approximation  of  the 


i>M"iM.il  lliiw  lii'lil  t>ii<l.  \  ll  III  iiiu.li  till''  wolilil  |>ti  i\  i<li' .1  trMi|iiii;it  \  liii<li><‘  lM'lv\<-4'n 

lli<‘  tliirrrnit  ^riil  roiiliKiiialioiis,  il  iiii^lit  ii«>(  Im*  a  Ioiik-Iciiii  Holiitioii.  To 


willi,  it  iin|tos«>s  an  adililional  Iturtlt'ii  upon  tlie  nsi'r.  StHoinlly,  using 
inloi polal ion  l«)  drvrUip  a  rorlangnlar  gii<l  and  llirii  inliM'polat ing  again  to  loratr 
largot  snrfa<  «*s  iniplirs  flir  polt’iilial  loss  of  a<'<  nra<  y  in  tin*  c^niisr  of  tin-  riMidning 
l»riM('ss.  Allliongli  \  Il'I'lU  is  inti'iidocl  as  a  *pialif at ivr  I«h)I  and  «an  I ln'i»‘f»)i(' 
lolnato  soiin  ina<  <  \na<  y.  it  is  <  Iraily  drslrald*'  to  limit  tliis  l«)ss  of  pns  ision  as 
mni'li  as  possildo. 


['ignre  12.  Snhvoinmc  Isolation  \  ia  Slab  llrinoval  Mctliod. 


\  Il’K.It  ran  now  inocoss  inrgidavly-si/»'i|  voxi  ls;  I  ilon't  Urlirvr 
n<»n-i«M  tangnlar  vo.\<‘ls  will  nspiin’ any  major  «  |iam>rs  lo  tin-  liilinrar  inloi polal ion 
or  volume  tra\‘<*rsal  routines.  However.  I  flou  t  see  an\  allernal  ives  to  storing  the 
.\.^  .Z  r»K)r<rmates  for  each  mesh  nofle.  If  this  provfs  nen-ssarv.  the  iiK-mory 
rerpiirerl  to  store  or  represent  the  data  file  will  l>e  rpiadriipled;  this  may  make 
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Figtin*  i;i.  Tyiiital  Arrtxlynainirisf ‘x  drills 

N'll’KH  imprarilral  for  iw  on  many  systnns.  Tlior«*fore.  I  wonM  Miiggrst  1hi» 
proliicni  l»r  <’N|»lonil  in  an  ♦•(fori  to  prwluii’  a  |>r<w(*«liiro  that  wonlil  allow  a 
griHTalizril  i|<‘n('i  iplion  of  Mu*  Iransforni  niininil  to  nia|>  grills  from 
non-ri'i  tangnlar  to  rts  t angular  fornis. 

jf.J.J  ,\oii~Sr»lnr  \  alufs.  AlVhongii  Ot)  )>i'r  ri*n1  of  tin*  stamlartl  flow  lioUl 
\arial»l«'s  ar<*  vrrlors  coinpoin'iits.  V  IPKR  ran  only  ilisplay  scalar  tlata.  Vector  ilala 
niaif  Ik’  «lispla\r<l:  liowri  rr,  siin  c  only  one  vi'clor  cuinponcnl  conM  l>c  ilisplayc»l  per 
iniagi*.  I  In*  value  of  sin  li  images  is  liighiv  ijnesi  i«»naltle.  I  lierefure.  I  re<  onnn<'n<l  ll»e 
prolilein  of  generating  nsi'fnl  displays  »»f  \eclor  data  lie  I'NaminejI.  One  possil»ilil\ 
wonlil  l>e  to  cominne  t  he  I  liree  \elocity  \e»  tors  ( t  I  I  - 1  ami  <lisplaying  the 
magnitude  of  the  resultant  vector  as  a  scalar  value. 


I  I  I  'll  I II  nil  1 1  r  I  h  I  I  nil  ion  (  li  »sc|\  ril  In  I  hr  |  >i  i  il  >lriil  nl  li  is|  ihi  \  iim 

\«‘rlnr  vhIii«‘s  is  I  In*  ins'«l  to  )lis|>li«y  ilitia  that  is  ^l^•|•iv»•^l  from  tin*  stamlard  (low  lirhl 
valiirs  of  inossiin*.  nn*rg_v.  ami  v«*lo«  ily  v«sl«irs.  In  or«l«*r  to  |»rovi«l«*  this  capahility. 
a  lanj^iiai^r  |)i<M«*ssor  must  In*  ih'\rh»|i<‘d  so  that  tin*  nsn-  may  doscriln*  how  tin* 
th’sircd  |>araim‘t<'i  is  driivod.  Mv  hs’ling  is  that  this  foatiin*  shoiihl  hr*  d»*\<*lo|>«*«l  as 
a  s«*|>arat«*  inoKiam  |•ath^*l'  than  as  an  int<*Kral  part  of  \  Il’KU:  ot h«*iwis(*.  \  IIM'.U 
may  l)«*«om«*  loo  ('mnlH*rsom<*  to  nsr. 

/.,*./  Siipi  iTnnii>iili  rs  mui  I'twi  ssnuj.  As  |)n*\ionsly  nu'iil  ion<*d. 

\  ll’I'.U  tHk<*s  t<K>  long  1.0  g<*m*rat<*  imag(*s  to  hr  nsrri  for  intrrar  tiv**  virwing. 
Ilow(*vrr,  it  is  poasihir  that  if  instalh*<l  on  a  snprrrompntrr  or  hrokrn  np  into 
srvrral  paralh'l  prof«*sst*s.  thr  timr  r<*((nirrd  to  complrtr  an  imago  rouid  possibly  h« 
r<*<lmr<l  to  soronds  instrad  of  honrs.  Uostnuturing  \  IPICH  for  parallrl  proc«*ssing 
shouhl  Im*  r<*lativr|y  aimph*:  individnal  srgmriits  of  I  ho  viowing  plane  [t  .g.>  a  fow 
scan  linos)  rouid  ho  assignoil  to  dilforont  prorossors.  When  ruiisherU  (ho  s<*gmonts 
wonhl  (»o  roimitod  to  form  tin*  final  imago.  'I’ho  pot«*nlial  for  improved  performam  r 
<h*(ini(«*ly  makes  this  area  worth  <*.xploring. 

Should  VIPKH's  »*N«*<  ntion  time  he  shortenerl  enough  to  permit  inl«*rarl iv** 
vi(*wing.  n<*w  user  interfaces  will  ims*<I  to  hr*  «|ev«*lope«l.  'rh«*sr*  interfaces  will  have 
to  allow  tin*  user  to  rotate  the  volume,  changr*  the  viewpoint,  ami  isolate  various 
suhr«*gions  of  tin*  volume.  Il  would  also  he  nsefnl  to  allow  the  user  to  seh*c( 
rlilferent  targr’t  vahn*s  or  ev«*n  riiange  the  type  of  <lata  h<*ing  displayed  (c.«/. .  from 
vehxily  l«>  pressnn*)  “on  the  fly". 

»  .\ iilomnlir  f'lnir  l  iihl  Inihm  Ihlnhon.  One  of  \  ll’l'df  s 

shortcomings  is  that  tin*  us«*r  must  ha\<*some  idea  ol  what  data  values  will  l>e  ol 
inl»*rr*sl.  If  tin*  usr'r  do«*sirt  aheadv  ha\«*some  idr*a  o(  what  (In*  (low  (i<*hl  should 
look  lik<*.  In*  must  r<*l\  on  r  hami*  to  find  sin  h  (ealnres  as  shook  waves  ami  vortices 
()m*  way  to  hnatr*  these  h'alures  is  to  pr<*prc»ress  the  data  to  (imi  minimum. 


III.)  \  iniiiiii  .  >1111 1  ,iM  I  \  ,1 1  Ill's.  \\  il  |i  till-.  I II  III!  1 11,1 1  lull .  I  III-  iisi'i  I  I  mill  I  riH  li'i  ,iii 

iisinn  vhIucs  ami  wiilr  nioiif^li  wimlows  to  sjiaii  all  of  I  In*  ilala.  Willi 

siirrcNsivo  iiiiaKi's,  i  rilii  al  How  (ioM  f»*at tin's  may  Iw  frirolisl  oiil.  1  lii.s  is  <  l«*aily  an 
awkwanl.  In  iilo  foirr  a|i|»roai  li  lo  an  im|iortanl  issin'.  It  wonhl  In*  Ix'llor  if  \ 
l  onM  anioiiiat ically  locatr  ami  rcmlor  arras  of  |)olrntial  inlrirst  lo  rngim'i'is.  I'liis 
fralnir  miglil  lakr  I  hr  form  of  having  tin'  nsi'i  sprrify  I  hr  |iaianirtri  lo  In' 
iiioiiilorril  ami  a  llirrshohl  for  I  hr  loc  al  graclirnl.  If  \  IIM'.K  ch'lcs  Utl  a  loc  al 
gradirnl  lhal  rxccs'di'cl  this  llirrshohl.  il  would  rc'iidri  I  In'  rrgion  ami  save*  ihr 
|>aramc'lc'r  value's  in  a  sr|)aralr  (ilr  for  |*osl-rrmlrring  analysis.  Snch  a  frainrr  could 
save'  a  gre-al  de'al  of  lime'  ami  frnst ration  for  llm  nsrr. 

^.J.6  Tniipnivl  /‘Vfl/uce.s.  Many  How  lirld  fratmrs  snch  as  shock  wave's 
and  vortiers  drvrlo|)  ami  change'  as  a  fnnctioii  of  limr.  It  would  t  herefore'  he  useful 
to  In'  ahir  to  s|n'cify  a  start  and  stop  time  along  with  the'  nnmlx'r  of  image's  to  he 
proelncrd  and  have  \'IPl'dl  aniomal ically  dc'Ve'leip  the  de'sireel  seepieme  of  image's, 
rids  capahilit y  would  obviously  re'epiirc'  the'  use'r  te>  supply  some  nmans  of 
a.ssoe  iat ing  a  spe'e  ific  time'  fe»r  e'ach  data  semiple'.  If  lids  coulel  lie  done  in 
eeuijiincl ieui  with  supe're'ompnt ing  ami  paralle’l  proce-ssing  (sen'  alxive),  the 
peissihilily  eif  auleunalic  ‘•men  ie'”  ge'iirralieni  ami  intrractive  image  vic'wing  might 
he'e'omr  very  re>al  indes'el. 

.f.J.7  Siibrolunn  holntion.  \  I PKR  allems  the  nse’r  to  isolate  data 
suhvolume’s  hy  explicitly  re’inoving  laye'is  from  the  voinmr  faee's  (which  is  much  like 
pe'e-liug  the'  layers  off  of  an  oide)u)  ami  rmhrelding  the'  vie'wing  plane  vvilhin  the 
data  volume'  (wide  h  is  like  making  one'  sliie-  in  an  emiou).  Ilowrve'i.  use-rs  will  ofle  ii 
waul  to  re'inov  e  se’e  lions  which  elei  not  lie*  paralh'l  to  I  In-  volume'  fac  e's.  Insle'ad.  I  he'v 
will  iie'c'd  I  he'  ahilil  v  lo  carve  oul  anv  ariiilrai  v  suhvolume'  I  he'v  so  eh'sire'.  l  ids 
e  apahilil  v  should  he'  prov  ide'd. 
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Mn.Niii*'  Hiowii  rrmilly  <»l»sri  v<*«|  tlial  llinr  ar«Mi‘l  vn  v  inaiiv  visiializal ion 
tools  lliat  run  on  woikstal  ions.  “Mach  scii'nlisl  is  inv<'sl  i^al  iiif;  sonic  lit  1 1<‘ as|i<'(  I  of 
some  fai -fct<hc<l  |>i'oli|cm,  so  they  are  not  i^oinK  t»>  liinl  general  tools’*  sh<‘  notes. 

"In  fact,  that  is  a  |>rohlem;  there  is  a  h'vel  wln'ie  tln'V  all  n<s*(l  similar  tools,  hnt 
then  past  that  level,  they  all  inssi  a  Rieat  (leans' of  spec  iaiizat  ion  (‘i'liJ’’)).  I  have 
ch'sif^nc'cl  \  Il’KU  to  o|)('rate  on  a  varic'ty  of  workstations  and  systems.  I  also 
developc'cl  \'IIM'.|{  in  a  general  mannc'r  so  that  spcs  ializc'd  applications  aiul  gc'iicMal 
c'ldiancc'ments  ronhl  he  easily  incorporatcnl.  Although  the  primary  ohjectiveof  this 
tlu'sis  was  to  produce  a  program  that  could  generate  images  of  a/'iodynamic  flow 
field  data,  volunu'tric  rendering  tc'chnUpies  lend  themsedves  to  a  wide  variety  of 
possible'  u.ses.  In  the  course  of  this  thesis,  I  have  rendered  images  of  ac'iodynamic 
(low  fields,  c'l  ror  levels  in  magnetic  fields,  and  atomic  energ,v  levels  hetwi'en  crystals 
without  modifying  \’II*KR.  'Phis  h'ads  me  to  conclude  that  VIPMH  cpialifi<»s  as  a 
“general  texir.  I  also  believe  that,  if  necessary.  X’IPKH  could  be  easily  modific'd  to 
ac  comodate*  spc*c  iaiizat  ion  in  various  (ic*|cls.  some*  (*f  which  arc*  di.se  usved  bc*|ow. 

i-i.l  liiidar  ('loss  Sirlioiis.  Radar  cross  s«'ctio\>s  arc*  an  important 
consideration  when  dc'signing  aircraft  anci  missiles.  espc*cially  at  a  time  where 
emphasis  on  “stc*alth  (<*chnology”  is  incrc*asing.  Since  radar  cross  sc'ction  data  is 
collectc*cl  in  a  dl)  mc*sh,  \  IIM’R  cocdcl  be  us<*cl  to  a.ss<'ss  radar  retunts  from  many 
look  anglc*s  as  wc*ll  as  diffc'ient  radio  frecpic’iicic's  and  power  settings. 

/  riinnl  I  '.nvf  lojK s.  W  het  her  em|'!ac  c'd  around  an  area  or  point 

targc’t.  r  <la:  cc)vc*ragc'  is  afrc*ctf'cl  by  the*  targc't’s  sinroniicling  tc-riain.  .\lso.  there* 
may  be*  zones  of  little  or  no  radar  covc'iage  bc'tvvc'c'n  installations,  bikc’wisc*.  anti-air 
clefc*nsc*s  (SAM  ami  ,\,A.\)  are  susceptible*  to  the*  same*  tc*irain  masking  and  site* 
separation  cc)nsiclc*rat ions.  I  hc*  air  combat  are  na  is  somc*\vhat  nnic|nc'  in  that  it  is 
t hrc*c'-climensional  instc*acl  of  the  t wo-dimc*nsion  rc*alm  as.soc  iatc*c|  with  ground  or 


ti.i\  .tl  -^lll  III  ('  1 1  iimImI  .  I  I II 'I  I'll  III',  it  ''I  It  III  III  III'  |ii  issil  ill'  III  Mil  II  Id  I  III'  V I  linn  11'  ;ii  uiiiiil 
ii  Kivi-ii  liHKrl  to  |)ro<liif  (>  hii  iiiitiK«'  Uial  iiU'iililioM  |>ot<’iil iaily  vulii«'i'al>l«‘  a.<<  \v«>ll  ax 
liil^li  lliii'al  I'ogioii.s.  \'II*KI{  <'o(il<l  ho  in  sikIi  a  smiario  ax  an  aid  to  nns.siun 
|)laniiiiiK  and  risk  a.ssrs.sm<'ut.  \  II’I’U  roidil  also  ho  iisod  as  an  iidolliRoino  Ux)l  if 
iisod  to  analy/o  ki'y  si{^natni'<‘  oloinonls  or  |>orform  s1r«’nglSi/w<*akiioss  ass<‘ssnK‘n(s 
»»f  woapons  .sysloiiis  siicli  as  oiioiny  fiftldors,  naval  surfaco  (oiiihalaiils.  or  tanks. 

{.  i.  i  luirlh  SrifiKts.  \'oluiiiolri<'  roiHloring  looliiii(|nos  havo  alroady  Imsmi 
appliod  (o  nioltsmilogical  data.  Si'voral  possihio  ii.sos  for  \  IPRR  may  also  Im*  found 
in  (ho  oarth  soionoos.  Oiio  aroa  of  rosoairli  (hat  could  make  use  of  VIl’KR  involves 
I  he  study  of  iso(h<‘rmir  layers  in  the  ocean.  yMthungh  generally  understood,  (here 
is  still  no  way  of  predicting  (heir  (hicknes.s  or  extent.  1'hese  layers  provide  excellent 
concealment  for  suhmariiies;  if  the  ocean  surrounding  a  surface  cumhatant  could  he 
samph'rl  in  .several  plac<*s  and  I  In*  resnllaiit  data  was  |>roc<*sse<l  l>y  N’lPlill,  a  31) 
image  of  the  isothermic  layers  in  the  imiiKHliate  region  could  l>e  produced;  this 
would  allow  naval  commanders  to  narrow  <lown  areas  where  enemy  suhmarim's 
coiihl  he  hi<liiig.  It  woidd  als<j  allow  deter  tion  (Hpiipment  o|)erators  to  rine'tune 
their  e(|uipmenl  to  lompeiisate  for  a  given  isothermic  conliguration. 

Allot h«’r  earth  science  application  for  Vll’KH  conhl  he  found  in  the  realm  of 
seismic  data.  Information  from  core  samples  ainl  seismic  sensors  could  he 
iiitegrat«'d  and  rendered  to  yield  a  31)  image  of  variotts  features.  Applications 
would  include  su<h  diverse  activities  as  earthquake  prediction,  petroleum 
e.xplorat  ioii.  and  assessing  the  elfe<  (  ivene.ss  of  hardening  terhni«iiKs  around  missile 
silos.  himk<'rs.  and  other  nndei ground  strucliires. 

/../  .Siinimarif 

Wendering  images  hase<l  on  \<»liimelri<  data  is  a  rela(i\ely  new  field  of 
einh'avor  in  computer  graphics.  For  this  (In'sis.  I  have  develop»'d  \  Il’KR.  a 


45 


v  |i<'  \  til  III  III 'I  I  i<  It 'll  I  If 'I  mu  |iitit;i,ini.  Mlt'i  )•ltM|l|l  mi^  i| )  iiii.ii^t-s  nl 
rin'oilyiiainit'  How  lirld.  Hloiiii<  niorgy  Irvols,  and  iiia|^n«*lii'  sriisor  data.  I  am 
(•oiiviiM  ivI  I  hat  V'll’KK  may  ho  iis»*«l  as  a  |M>wrrfid  rnRiiHHritiR  tool.  Tlu'  romnu’iits 
nM  oivrd  from  sovoral  irotoiilial  ns<'rs  s<'«Mn  to  «<iiilirm  this  hrlh'f  and  justify  more 
dovolo|im«‘nl  «*frort  in  this  arra.  AlllioiiRh  vohim«*tii«  n-rnh'iiiiR  rannol  possihly 
r('l>la< o  ot hor  imaR<'ry  t«M»|s  sm  li  as  int«‘ra<  tivr  partirh*  Irarrs.  I  hs'l  that  it  is 
d<'stin«'<l  to  tak(‘  its  place  as  an  augmentation  to  the  array  of  aids  availahh'  to 
aerodynamicists  in  tin'ir  pursuit  of  safe,  efiicient  aircraft  design.  \N  Idle  (h’velopiug 
N'IPKH  ainl  analy/ing  it.s  applicahility  to  aeroilynamic  engineering.  I  have  InH'ome 
aware  of  the  range  of  possible  us<*s  for  programs  such  as  \’IPKIt..lt  is  my  belief 
that  as  each  area  is  explored,  more  use's  will  be  runt  iniionsly  nnrovereil.  Therefore. 
I  be'lieve  that  VIPKR  shoidd  be  enhanced  ainl  its  applicability  to  ae'rodynamir  (low 
fie'Ids  as  well  as  other  areas  of  research  and  developtiM'iit  should  be  |>nrsuetl. 


Apix'mlix  A.  Ust-r's  Gault 


A.  I  Intimiurlion 

riu'  VoluiiuMiic  I’rogiam  for  Kogiucoriug  U«'sfan  li  (V'll’KH)  is  a 

ftriK'ial  |>iir|)osf*  gra|ilii('s  program  (l(*sigii<Hl  lo  pioHm  r  tlir«H'-(limi*iisioii  (31)) 
imagos  Uasrfl  on  m«'s|i«>s  of  volimirt  rir  sralar  data.  \  IPKK  is  srini-iiilorarUvr  in 
that  it  allows  tin'  iis<'r  to  interact ivoly  s«*t  the  parameters  that  control  tin*  image  to 
In'  pro<lnee(l. 

This  User's  (!ui<le  will  outline  how  to  install  VIPEK  on  a  typical  computer 
system  and  how  to  (h'termiin'  the  parameter  values  you  will  need  to  prorluce  the 
image  i/oii  want.  Finally,  some  hints  for  .speeding  up  image  production  will  he 
providtxl. 

Although  N'll’ER  was  develop«'d  on  a  Uni.x  operating  system.  I  have  tried  to 
k«s'p  the  program  as  system- independent  as  possible.  Therefore.  1  Ix'lieve  that 
\  ll’EH  should  run  on  any  system  that  suppmts  the  programming  language  and 
has  sulficient  memory  to  handle  large  data  Hh-s.  Likewise,  this  Us«*r's  (iuide  is 
1/iased  towards  Uni.x  syst«'ms;  however,  outsi«le  of  some  minor  installation  details, 
there  should  be  no  real  complications  involv«sl  with  using  VIPEH  on  non-Unix 
systems. 

.I.J  Installing  \'IPI:'H 

\  ll’EH  actually  retpiires  two  librarh's  in  order  to  operate.  Phe  first  library 
is  c»)mpos<'d  of  tin’  fih’s  that  make  up  tin*  Utah.rle  imag«'  foiiuat  loutines.  I'hesc' 
routines  are  a  public  domain  set  of  files  that  <'nable  \  IPEH  to  pio»luie  mor«' 
e((ici<'nt  image  files  as  an  output.  I  he  se<<)nd  \  ll’EH  liluarv  is  compris<'<l  of  tli*' 
files  that  produce  \  ll’E.H  itself. 


I ,  /  I  lull  I  li  I  I  hi  II I  II  1  .il  ill'  I  IinI  s  I  In-  files  I  liitl  liiftki’  11 1 1  I  ln‘  I  ( ;ili .  I  Ir 
lil»i'«fy.  SoMH'  of  I lii'se  files  in«v  in'*'*!  lo  Iw  iniMlifieil  to  Dpernli*  on  your  coinpiiler 
system.  Inst riiel ions  for  iixxliiieMl  ions  may  Im*  foiiiHl  within  the  files  aikI  are 
heyoinl  tin*  scope  of  (his  user's  f^ni«le.  Once  tin's**  fib's  have  Ins'ii  installeil,  yon  will 
IKS'*!  to  eompih'  ainl  link  tln'in  to  form  (In'  final  library  iinnlnh'  librie.fi.  The 
simph'st  way  to  <lo  this  is  to  n'liann'  rleinakefile  (*>  makefile  and  (hen  «'n(er  the 
eommaiid  iiiaki .  Doinj*  this  will  <  ans«'  makefile  to  antomal  i«  ally  pnnhire  librie.a 
for  y*)n  {.\oh:  This  applit's  to  I'nix  systems;  if  yon  are  a  non-l'nix  n.s<'r.  makefile 
may  not  work  for  you.  Consult  your  .system  manager  or  dormnentation  to 
d«'termine  the  compile/link  op«'rations  you  will  ne<^d  to  perform  for  this  st<*p.). 
Once  librie.a  has  In'en  prepare*!,  you  are  ready  to  install  VII’KH. 


bst  ring.*’ 

rb’^*’tr*»w.*- 

rb’niak*’file 

sv-pntraw.c.fi 

buibimap.c 

rb’^etrow.c.b 

Kuiisv.c 

sv.pntrow.c 

dither. c 

rle^etrow.c.f 

scanargs.r 

svn)^lobal.r 

rle.getroni.c 

rb’_put<om.c 

svfl>.h 

svn)..gb>bal.h 

rle.g*’traw.c 

rb’jaw^lc.c 

sv.putraw.r 

XindRunsv.h 

rle.gf't  raw.h 

rlej-ow-alr.r 

sv.piitraw.r.ba 

Table  I.  ('tail. lie  Library  Fil**s 

A.J.J  \  ll’l\H  l.ibivry  The  (ih's  y*)u  will  nrs^d  for  N'll’KH  are  shown  in 
d’able  2.  Once  tin’s*'  fib’s  are  in  pla<e  <»n  your  syst*’m  ainl  you  ha\e  prodiir*’*!  a  copy 
of  librie.a  (s*’e  above),  you  ar*’  ready  to  Imibl  the  V  Il’KR  library.  To  do  this, 
rename  vipermakefile  to  makefile  and  enter  the  eommand  makf  riinv  (*> 
romplete  the  \  Il’KH-buibling  s*’*pn’nre  (again,  non-l'nix  us«'rs  may  ins’tl  to  follow 
a  diffi’ix’iit  pr*)<'*’*lnr*’. ).  ^on  should  now  hav*'  ail  \'on  n<’<'*l  to  produ*'*’  'ID  image’s  of 
your  v*)lnm*’t  i  i*'  *lal  a. 

A.  i  ('Sing  \  I ri:H 

('sing  \  11’1‘dl  is  a  fairly  straight  forwanl  proe<’ss  that  n’»|nires  a  lib’ of 
properly  fiuniatted  scalar  data  and  a  f*’W  u.ser-supplied  parameter  valiif’s. 


■ 


F 


r. 


■ 
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I  nlnl  <  misi  .<  s  li 

•'Xlrriis.li 

iiifrriMtIalr.r  liavrrHal.f  vl|M'r.li 

inapvHl.r  v«‘x(<‘riiH.|i  vi|><‘t'inakf‘fil(' 

mat  rix.o|)s.(- 

ral»li*  2.  \  UImhi  v  l  il.'s 


A.  i.l  Ihtin  hhnnnl  I'iftmo  1  I  slums  a  small  lilr  of  ilala  in  tin-  \  IIM’.H 
format.  1  ho  first  lino  is  simply  a  tomarks  lino.  'Phis  is  a  lino  of  np  to  SU  chaia(t«Ms 
ami  is  rosorvod  for  any  lalnds  or  portinont  information  that  idontilios  tho  fd**.  Tho 
contonts  of  t  his  line  are  immaterial  to  VIPER:  however,  this  line  must  be 
present)  even  if  it  is  blsnk.  1  he  secoiul  lin<‘  contains  thrrx*  inte^i'rs  that 
imlicato  the  (limotisions  (in  terms  of  data  nodes)  along  the  X.  ) .  and  Z  axinr  of  the 
data  Volume  to  he  ronderetl.  'Phe  remaining  lines  are  comprised  of  the  A’.  Z 
coordinates  of  each  data  point  along  with  the  scalar  <lat  K  value  at  that  point  {Xoir: 
Putting  OIK'  s<*t  of  coordinates  and  value  |)er  line  makes  tho  file  easier  to  read; 
liowever,  the  file  will  rerptire  more  memory  space  if  you  do  t  his.  For  larger  fih*s.  it 
may  Ix'  hetter  to  r«'cord  iuok'  Ilian  one  data  point  p«'r  lim’.  \  IPER  ran  read  eillu'i 
format.).  Hememher.  VIPER  assumes  that  you  are  using  a  n d a ngulai  grid 
of  data  points!  ,\ny  other  type  of  grid  will  he  distorted  and  the  final  image  will 
he  dillicult  (if  not  impossihle)  to  interpret. 


Sample  VIPER  Data  File  -  Comment  Line 
2  2  2 


-2.U 

-2.0 

0.0 

-l8.tKH) 

-2.t) 

-2.0 

1.0 

12.000 

-2.0 

-1.0 

1.0 

C.OtKt 

-2.0 

-1.0 

2.0 

S.IMUl 

-1.0 

-2.0 

1.0 

(i.OOO 

-1.0 

-2.0 

2.0 

o.ttoo 

-1.0 

-1.0 

0.0 

0.000 

-1.0 

-1.0 

2.0 

:{.ooo 

Figme  II.  Sample  V  IIM'.H  Data  File 
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1.  J  I'iiniiiii  U  r  Silhmi  \  IIM'.lf  will  \<im  In  m-I  scvciiil  \i<'winj> 

imraiiK'liTs.  In  *mlrr  fur  you  to  In*  al>lr  to  iiiakr  inforiii«‘tl  valiu'  (wlrrtioiis,  vve  will 
now  r»*VM'\v  «'h<  Ii  iiaranK'tri'  ainl  liow  (o  «lH<*rniin(‘  tin*  vhIik's  yon  will  n<H'«l. 

A.  i.J.I  ( 'onininml  Lint  l^'nlrits.  'I'lit' lirst  |>arHinHiM's  \'on  will  sH 
aro  niil»f’<l«l(‘(|  in  tin*  roinniaiul  lint'  iimhI  to  invuko  VIPICIi.  I'lu'  eoniinand  linr 
lornial  is 

viper  -option  input  output 

'I'Ik'  -option  allows  you  to  locate  the  |K>aition  of  the  image.  Using  -c  will 
center  the  image  on  your  screen  while  -k  will  locate  the  image  in  the  lower  left 
corner.  Normally,  yon  will  probably  want  to  center  the  image;  the  -k  option  was 
(l<'v«'lope<l  for  s|)ecial  situations  (such  as  producing  vi(leota|H*s  of  images). 

The  input  parameter  is  the  name  of  the  data  file  you  wish  to  have 
remh'red:  output  is  the  name  you  assign  to  the  final  image  file  {^'otr:  This  file 
name  must  lia\«‘ a  .rie  «'.\l<‘nsion.). 

If  your  command  line  doesn't  have  the  proper  formal,  VIPBR  will  reject  it 
and  terminati*.  Otherwise,  VIPER  will  reail  your  data  file  and  then  enter  the 
interactive  parameter-setting  phase. 

A.  LJ.J  Ffttihack.  The  first  value  yon  will  be  asked  for  will 
determine  whether  or  not  VIPER  produces  feetiback  for  you  during  image 
generat  ion.  This  feedbai  k  takes  the  form  of  echoing  back  some  of  your  parameter 
selec  tions  as  well  as  printing  an  advisory  notice  afti'r  eac  h  siaidim'  in  yom  image' 
has  Ixvn  completed.  I'liis  information  is  useful  if  \«>u  want  to  monitor  th*'  progrc'ss 
of  your  image*  as  it  progresses  as  we*ll  as  the*  first  fe'w  time's  you  use'  \  IPh.R.  Once' 
you  get  list'd  to  rimiiing  \  IPER  eir  want  to  put  \'IPER  iiite)  the  baekgroimd  mexle. 
ye)u  will  probably  want  te>  elisabh'  the'  le'e'elback  option,  lo  e'uabh'  fes'dbae  k.  enter  1; 
to  elisabh*  h'e’dback,  ente'r  0. 
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I  i'.  i'  Sr)ilin</.  Allri  \  IIM’.U  liiis  hm<I  vuiii  il;i(.t  lilr.  il  will  Irl 
yiHi  know  llir  tiiiHl  voliiiiH'  tlitiinisioiiH  aii<l  ask  yoti  foi'  a  si  aliiif^  lactoi'.  I  Ills  is 
iipport  unity  to  rltangr  tin*  pliysiral  (liiuiuiKionM  of  your  voluiiic  along  any  or  all  of 
t  in*  voinuM'  s  axi*s.  Ho\vov«‘r,  to  tlrttuuiim*  tlu*  proper  vaitio  will  r»*<pnr«*  you  to  set 
several  paraiiu'ters.  In  or«ler  to  work  out  these  values,  we  will  now  look  at  the 
viewing  g«'oni«'try  that  VII’MH  uses  (s«h‘  Figure  l.'i).  N’lPKIl  assumes  that  you  are 
at  the  axes  origin  and  looking  straight  <lown  the  Z  axis.  I'he  viewing  plane  and 
data  voluiiK*  ar«*  initially  eeutei«‘d  on  tin*  ^ axis.  Vll’KIl  will  allow  you  to  set  the 
distanees  from  your  eye  to  the  viewing  plane  as  well  as  the  centiT  of  the  tiata 
volume.  This  capahility  allows  you  to  use  VIHKR  in  murh  the  same  way  as  j-oii 
wouhl  use  a  ranu'ra.  Moving  the  viewing  plane  is  inueh  like  adjusting  a  camera’s 
focal  length  while  moving  the  volume  corresponds  to  moving  towards  or  away  from 
an  oljje<  t  in  order  to  adjust  the  total  area  It  will  occupy  in  a  photograph.  However, 
V’IPKR  gives  you  one  more  option:  placing  the  viewing  plane  insidt  the  volume 
allows  you  to  cut  away  port  ions  of  the  volume:  whatever  is  between  your  eye  and 
the  viewing  plane  will  be  omitt<'d  from  the  niial  imag«*. 

Referring  to  Figur<'  10.  we  see  how  to  determine  the  necessary  parameter 
s<*ttings  by  using  similar  triangles.  We  already  know  the  si/e  of  the  riata  volume,  d. 
If  we  then  choose  a  value  for  the  distancr*  from  our  eye  to  t  he  viewing  plane  (n), 
the  desire«l  extent  of  the  volume  in  the  final  image  (6).  and  the  distance  from  our 
eye  to  the  anfrroi  the  <lata  volume  r.  we  may  ea-sily  fitid  the  scaling  factor 
refpiireil  to  [Moduce  this  image  by  the  following  ecpiation: 

ScaU  Forfar  = 

'I  X't 

Once  you  have  determineil  the  scaling  fa«  lor.  \  t»u  shonhl  eiilcM  il  into 
\  IPFR.  If  you  do  not  «’ntei  tin'  same  vahn*  for  all  three  axes,  tin'  data  volume  will 
!»<■  “strc'lched"  from  its  original  proptutious.  .Vuother  thing  to  k«s'p  in  mind  is  that 
if  you  <lon  I  want  to  change'  tin'  volunn'  s  si/<'.  you  should  <'iil*'r  a  value'  of  l.lt. 


51 


Figure  l(>.  CaUulHtiiig  Oljject  Positions  and  Sealing  Faetors 

Mauiimlaliug  th(  Data  S'olamr.  VIPER  will  next  give  you 
I  lie  <)|>|)oi  I  unity  to  i<H>rien(  the  data  volume.  Ilie  lirsf  set  of  parameters  will  allow 
you  ti)  rotai»*  the  \'olum<'  about  any  or  all  of  its  tlirei*  axes.  Those  valiu*s  arc 
iiu’asurerl  in  degrei's  ami  in<  reas<'  in  a  anti-elo«  kwisr  diror  tion  (as  you  look  along 
the  axes).  'I'lie  second  set  of  parameters  permit  you  to  “slitle"  the  volume  along  the 
thrrs*  axes.  I’o  move  h'ft.  <lown  or  towards  you.  you  will  m*ed  to  enter  a  negative 
valm-:  likewise.  p«isi(ive  values  will  move  the  volume  right,  up.  and  away  from  you. 

.  l..y.  J.  7  Changing  Iht  \  it  wing  IVant  Dimmsions.  The  viewing 
plan<>  is  measured  in  terms  of  picture  «'lem«’nts  (/>/.ri/.s).  The  default  Inu'iT’.ontal  and 
vertical  dimensions  for  the  viewing  plane  are  1021  x  TiiS.  respect ividy.  I'o  maintain 
tlu'se  values,  enter  a  0  wlnsi  \  IPF.U  asks  you  if  \on  want  to  change  the  scrcsni  si/e. 
1'o  <  hange  these  dimensions,  enti'r  a  1  followed  l>y  the  new  hori/ontal  and  vertical 
dimensions  yon  want  (.\nli :  If  the  \alm*s  you  enter  exceeil  the  default  values  or  are 
less  than  two.  VIPFR  will  revert  to  the  default  scre<*n  size.). 
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i.  I.  '  I  iijlils.  I  lie  iicxl  M-l  ol  v\i||  i  |i 'I  < 'i  IllilX'  llu' 

liKhtiiiK  williiii  your  somm*.  I'lu*  first  valin'  r*‘|m‘snils  llu*  aiiiliinit  11^1)1  prosr'iit. 

I  Ills  is  till*  aiiioiiiil  of  IIkIiI  llial  \voiil«l  Ih-  availalili*  if  iioiir  of  your  liglit  soiirros 
won*  tiiriusi  on.  A  valm*  of  0.1)  corrrs|)uii(ls  to  a  totally  darkciird  room  wliilr 
I'lilrriiig  1.0  would  iiidicatr  a  liriglit..  siiidit  sfoiu*.  Noxt.  \'1PKH  will  ask  you  liow 
many  light  soiircrs  you  waul  to  install,  ^'ou  must  liavr  at  loasl  olio  light  hut  no 
iiioro  than  fivo  lights.  I'or  oarli  light  source,  you  will  ii«‘ed  to  supply  X’ll’KH  with 
the  tliris*  r(M>rdiuates  for  the  light's  location  (Heniemlier.  your  eye  is  at  the  origin 
with  a  locat  ion  of  (O.O.Oj.  A  light  100  units  rlirectly  above  your  eye  would  have 
coordinates  [0,0,100].).  When  all  of  your  lights  have  Ixvii  activaterl,  you  will  iiecti 
to  toll  VIPElf  what  it  is  supposeii  to  l>e  looking  for. 

Selling  Ihe  Taegt!  \alutA.  \'IPICI{  will  try  to  find  surfaces 
whoso  \aluos  corrospond  to  the  targets  you  s«*t  for  it.  Non  may  specify  from  one  to 
five  target  values  to  lie  sought.  After  you  tell  X’ll’ER  how  urauy  values  are  to  l>e 
found,  you  will  need  to  provid«*  a  set  of  values  for  «‘ach  target  you  select,  riie  first 
parameter  is  the  value  of  the  surface  itself,  ^'ou  may  then  8|>ecify  a  “window"  of 
values  that  are  “close  enough"  to  the  target  value  {(.g..  if  you  specify  a  target 
value  of  5t).(J  and  an  allowalih*  variance  if  o.t).  the  window  of  values  will  range  from 
15.0  to  5.5.0). 

I  he  lu'xl  |)aram<’ters  will  determine  the  surface  characteristics  of  color  ami 
opacity.  \'11’ER  produces  color  l)y  additiv<*ly  blending  red,  grer’u.  and  blue  (IKJB) 
to  yh’ld  the  <lesire<l  color  ( A'o/f;  ndflilirtlg  means  that  lower  color  valm*8  are 
•  larker;  thus,  an  ItflM  valiu’  of  000  pnalmes  black:  l))l|  is  blue.  I  I  I  is  while.  el<  ). 
kor  each  surfac**.  \ou  will  need  to  pro\i«h'  an  IKIM  triplet  to  spec  ify  th<’  color  of  the 
surface.  Be  careful  not  to  give  the  same'  IU!M  value*  to  more*  than  one  taigc't:  l>oth 
will  be  givc'u  the  same  color  and  you  won't  be  able*  to  tell  the  two  apart.  I'inally. 
you  will  provide  an  opacity  coefficient  for  tin*  surface.  I  his  will  ch’tc'iiniue  how 
transparent  the  surface  will  be.  A  value  of  0.0  will  be  comph'tely  transparent  and 
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Will!  I  -'liiiw  i||i  III  llic  lllliiKi';  .III  (i|i.i(ilv  III  III  v\lll  Ik*  i'iiiiiiiIi'Ii-Iv  ii|i;ii|I|i'  ;iiii|  I  IiiI'n 

Hii\  iiii<l<‘rlviiiK  stiiTa«'i‘s.  Any  vhIik'  Im-(\v«'<*ii  iIk'm*  two  limits  will  |>roviil<'  a 
part iaily-l  rHii.s|>ar<‘iU  surfftre. 

A.J.J.IO  Surfort  lit Jltdirilif.  Thoir  is  one  linal  surface  parameter 
that  is  applie<l  uniformly  to  all  surfaces.  'I’lie  surface  redeclivily  of  the  surfaces 
rauj^es  from  0.1)  for  a  very  <lull  surface  to  1.0  for  a  very  shiny  surface,  ^'ou  will  he 
al)lc  to  (lelcM'iniue  the  glossiness  of  the  surface  hy  entering  a  single  value  within 
these  two  homuls. 

A.  i.'J.H  Image  Maitt.  Tlie  linal  parameter  yon  will  l>e  prompted 
for  is  the  IKIU  triplet  for  the  image  matte.  This  is  a  default  color  that  will  be 
applit'd  to  any  part  of  the  viewing  plane  that  is  not  occnpunl  hy  the  data  volume. 
'I'he  imag*'  matte  is  useful  in  a  number  of  ways.  To  l>egin  with,  the  matte  will  allow 
yon  to  s<'e  where  the  data  volume  is.  1'his  not  only  helps  you  to  fine  tune  your 
vi«’wiiig  gtsnnetry  in  subsetpu'iit  images,  but  also  helps  you  to  orient  rendered 
sm  la<  es  in  n-lation  to  tln‘  volume.  If  you  <loirt  want  an  image  matte,  simply  enb'r 
the  IKiM  triplet  for  black  (0.0.0). 

()nc«*  yon  have  entered  this  sei|nence  of  parameter  values.  VIPER  will  enter 
its  Hiitonomons  mo<le  aiifl  proce<*d  with  generating  the  image  yon  have  specified. 
When  VIPER  has  comph'ted  your  imag«‘.  it  will  ask  you  if  you  want  to  generate 
anotlier  image.  I'.ntering  0  will  cause  VIPER  to  terminate.  On  the  other  hami,  if 
Von  enter  1.  VIPER  will  retain  the  data  volume  and  scaling  factors  from  the 
original  im.ige  and  reenter  I  lie  parameter-setting  phnse.  'I'his  cycle  may  be 
rejieated  for  as  many  image's  as  yon  eh'sin-.  I  lms,  yon  r  an  save'  a  great  di’al  of  limr' 
by  not  having  to  rr'peaf  the  riata-rr’ading  portion  of  I  lie  program. 
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I  /'  Sjutilniii  /  ft  Inniiii  ( ii  iii  nil inii 

i'h«'i«*  HIT  ii  iiinnli«'i'  of  w«ys  to  slioiim  the  tiim*  r*Mniiietl  to  generate  an 
image  with  \  Il’KK.  Most  of  them  entail  s«»ine  sort  of  tra<hs)(r  hetwtvn  tlie  lime 
re<|nire<|  ainl  the  detail  and  «|nalil\  of  the  final  imagi'.  I'his  stslion  will  hrielly 
outline  stuue  t)f  these  t»'ehni«iu»‘s  and  any  jM>ssil)le  tradindfs  that  I  havt'  l>«'en  aUle  to 
ideni  ifv. 

I'he  lirsi  t('('hni(|ue  involves  vniting  tiown  the  size  of  the  viewing  plane  a.s 
muvh  as  iM)ssil)h'.  Sinre  the  titne  re<|nir«Hl  to  generate  an  image  is  direr! ly 
proportional  to  the  mimher  of  pixels  that  comprise  the  image,  rutting  *lown  on  the 
niimher  of  pixels  will  cut  <lown  on  the  lime  rerpiiritl.  I'his  is  a  fairly  good 
terhni<|ne  in  that  it  do<*sn't  necessarily  d«'grade  the  final  image.  Instead,  it  re<piir«*s 
you  to  carefully  size  your  image  of  the  volume  and  trim  the  viewing  plane 
accor<lingly.  Naturally,  if  you  use  a  smaller  image  size  for  your  volume,  you  will  l>e 
al)le  to  trim  the  viewing  plane  even  more.  Also,  smaller  image  sizes  will  prorluce  a 
fast<’r  run.  even  if  you  don't  trim  the  screrui. 

.Another  way  to  shorten  the  execution  time  is  to  trim  i*way  uiiiKHessary 
portions  of  the  data  volume.  As  previously  menlioiusl,  this  can  he  accompli,shetl 
«'ilher  hy  emhedding  the  viewing  plane  within  the  data  volume  or  explicitly 
defining  slabs  of  the  vohmie  that  are  to  he  trimnusi  off.  This  spt'eds  up  the 
reiulering  process  hy  pres<-nting  h’ss  <lata  to  he  rendi‘re»l.  ’ton  should  In*  v<'ry 
laieful.  Iiowever.  when  using  this  t<>chni(|ue.  Indiscriminate  trimming  may  cause 
\  ll’F.H  to  omit  important  features  within  the  volume. 

losing  very  dull.  opa<pie  surfaces  will  pioduc*-  a  sluulei  ruu  time.  However, 
ojvarpie  surfaces  will  mask  any  uiulerlying  surfaces.  Howi-w'i.  this  teclmivpie  is 
useful  for  ‘‘tiuick  looks"  and  situations  where  you  know  there  will  Im*  no  uuderlviiig 
surfaces. 

A  final  tevlinicpie  1  have  ns<'d  to  sjns-d  up  the  image  time  is  to  place  the 
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\  |i'\\  ini'  |)|.iiii'  \i'l  \  I  hisr  Id  I  lie  «>\c  |M>s|tidn.  I  Ills  w  ill  i  .iiisc  t  lie  i  ;i\  s  iisctl  |d| 
iiUHRi'  ii'iiih'iiiiR  It)  Im*  iluirt'  tlivt'igt'iil .  .\s  lliis  tlivfim'iui'  iiicrrast's,  tli«'  liinr 
r(‘<|iiiri‘<l  to  protliirr  an  itimgt'  •U'froa.st'M.  Al  ll»*  santr  tiino.  inoro  divorgt'nt  rays 


■ 


le 

■ 
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|>ro<lii<'('  ev«>r-iii(  roa.sing  gaps  wIumo  no  clala  is  iiu  lndiHl  in  tlie  final  iniagi*. 
SigniliranI  foalnri’s  (sncIi  as  vorlio's  and  shock  waves)  conid  lx*  conci'alctl  in  tin'sc 
gaps.  Once  again.  I  his  t(*chni(pie  is  more  sniletl  for  “tpiick  look"  iniagt's  or  images 
of  volumes  I  hat  contain  only  coarse  features. 

As  I  previously  mentionetl.  tht'se  teciini<|nes  will  sptssl  np  t  he  rate  of  image 
e.xecntion;  however,  you  must  he  careful  and  weigh  the  coiiwxiuences  involvetl 
In'fore  proceeding  with  these  or  any  other  shortcuts. 

.  I .  •  j  Sail) pit  Sf  .s.s/oii 

The  following  is  a  transcript  from  an  actual  V’IPKR  session.  The  holdfaced 
entries  are  VIPER  prompts;  user  inputs  are  shown  in  standard  i)rint. 

viper  -c  error. vol- 

Do  you  want  feedback  (0  =  NO,  1  =  YES)  =>U 
Polhenius  Error  Field 

Please  enter  name  for  image  file  ^sam|>lel.rle 

The  Volume  Being  Rendered  Has  XYZ  Voxel  Dimensions  of  5  5  3. 
and  Goes  From  -2.00  to  2.00  Along  the  X  Axis, 

-2.00  to  2.00  Along  the  Y  Axis, 
and  -1.00  to  1.00  Along  the  Z  Axis. 

Enter  X,  Y,  &  Z  Scaling  Factors  ^I  'l  |.')  I  'l 

Want  to  Trim  Away  Part  of  the  Volume?  (No:0  Yes;l)  =>  I 

Trim  Off  the  LEFT  Side?  (No:  0  Yes:  1)  ^  Ii 

Trim  Off  the  RIGHT  Side?  (No:  0  Yes:  1)  =>() 

Trim  Off  the  TOP?  (No:  0  Yes:  1)  -  ii 
Trim  Off  the  BOTTOM?  (No:  0  Yes:  1)  =t^() 
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'IVim  OfT  file  FHONT?  (No:  0  ^'es:  1) 

How  Much?  =>'■). 7 

Trim  Off  the  BACK?  (No:  0  Yes:  1)  ^0 

How  Far  From  Your  Eye  to  CENTER  OF  VOLUME  =>!  '■>()(). 

How  Far  From  Your  Eye  to  the  Viewing  Plane  =>IUIH). 

Enter  Angle  of  Rotatio::  About  X  Axis  (Degrees)  =i>0. 

Enter  Angle  of  Rotation  About  Y  Axis  (Degrees) 

Enter  Angle  of  Rotation  About  Z  Axis  (Degrees)  =>(). 

Enter  X,  Y,  and  Z  Iranslation  Distances  =»^U.  U.  0. 

Do  you  wai.»  to  change  the  screen  size  (0  =  NO,  1  =  YES)  =^l 
Enter  new  screen  dimensions  (x  y)  =>3r)0  3r)0 
Enter  Ambient  Light  Level  (0.0  — ►  1.0) 

How  Many  Light  Sources  Do  You  Want  ^2 
Enter  X,  Y,  Z  Location  for  Light  #  1  =>35.  50.  200. 

Enter  X,  Y,  Z  Location  for  Light  #  2  =^-35.  >00.  800. 

Enter  number  of  target  values  to  be  found  =>2 
Enter  target  value  #  1  =>0. 

Enter  allowable  target  value  variance  =>.5 

Enter  RGB  values  to  be  associated  with  0.00  =>  .5  0. 

Enter  Opacity  for  0.00  (Clear:0.0  — »  Opaque:  1.0)  =>.8 
Enter  target  value  #  2  =>2. 

Enter  allowable  target  value  variance  =>0.75 

Enter  RGB  values  to  be  associated  with  2.00  =>0.  0.  I . 

Enter  Opacity  for  2.00  (Clear:0.0  — ►  Opaque:!. 0)  =:•  I . 

Enter  Surface  Glossiness  (Dull:  0.0  — *  Shiny:  1.0)  =>  .!>5 
Enter  RGB  Values  for  Border  Matte  =>  1.  0.  0. 

Saving  File  as  sample!. rie 

One  more  time?  (  0  =;■  No,  1  =>  Yes  ):  0 
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1  .  ft  (  /ll.stflM 

(ii  (Icvflujfing  \  II’KU.  I  liAvr  tiiotl  to  pitxliu  c'  a  roiivoiiifiit .  rra.stMial)ly 
fftiririit  program  that  will  r*'inl«'r  U)  iniagt's  of  voliiiiiHric  data.  I  havo  tri«'«l  to 
iiu'orporatc  as  many  tiaofiil  foatmt's  as  t  t  ouUI.  If  you  can  think  (»f  any  tit  her 
capaliilit ies  that  would  he  usefid  (or  implement  any  of  them  yourself)  or  einountf'r 
prolilems  that  1  havt*  missed,  please  notify  the  computer  graphics  faculty  her<‘  at 
the  Df'partment  of  Klectrical  KngiiU'ering.  Air  Ktiire  |nstitut(>  of  recimology, 

\\  right- Patterstiu  Air  I'tirce  Uase.  Ohio. 
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(  apfHin  David  .lamrH  Bridges  was  Immi 

He  gifiditaled  from  Marion  A.  P<'f<'rs()ii  High  Selund  in 
After  marrviiijc^ 


birlli  of daiigliler, 
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ABSTHACT 

The  heliavior  of  fluids  has  been  shidied  for  years.  This  behavior  has  been  niodeled 
by  the  Navier- Stokes  e<|nalioiis  since  ihe  iiiid-iiineteeni h  century;  however,  tlii'se 
e((n!itions  are  so  roni|dex  that  their  use  by  engineers  was  impractical  until  the  advent 
of  th<‘  modern  computer.  The  same  computers  that  made  these  equations  iisefni 
A'4Screa(ed;a'uew  t>roblem>dHta  saturation.  The  solution  to  (he  Navier-Stokes  e<|ua(ions-, 
which  repri'sc'iit  an  aerodynamic  flow  field'wT^tuiiprised  of  (hoii.sands  rf-not  niillion^^f 
numbers.  This  much  datSiii-nkes  it  virtually  impns.sible  to  conceptualize  the  situation 
within  the  flow  fiehl.  To  alleviate  this  problem,  various  computer  graphics  t^'chniipies 
have  be<'n  nse<|  to  provide  images  of  flow  fields.  Some  of  these  techniques  rely  on 
approximating  isometric  surfaces  of  interest  wifli  polygons;  others  only  provide  an 
image  of  a  two-dinien.sional  “slice"  of  the  flow  field. 

.This  thesis  projeclMpplies  a  relatively  new  graphics  techniipie  known  as  “volumetric 
rendering"  to  generat*y  three-dimensional  iniage^of  aerodynamic  flow  fields.  -Vol¬ 
umetric  rendering  has  shown  promising  n-sults  in  several  applications  of  scientific 
visualization.  One  ^dvantageiof  this  t.echiii<|ue^*}hat  it  doesn't  rely  on  geometric 
primitives  to  approximate  a  surface,  ^AlSo,  no  a  prion  knowledge  of  the  flow  field's 
form  is  required  to  produce  the  final  ttnoge>  -  ■  ' 

To  as.sf'ss  th<-  utility  of  volumetric  rendering  in  the  field  of  computational  fluid  dy¬ 
namics  (CFO),  a  program  known  as  VIJ’ER  is  developed.  This  thesis  outlines  the 
requirements  aud  specific  algorithms  for  this  program.  Images  of  flow  fields  are 
prf'sented  and  di.scuissed  along  with  program  enhancements.  “Although  VU’ER  is  de-  ' 
veloped  for  CFD  purposes,  its  design  permits  its  use  in  4iAerent  applications, sc'veral 
of  which  ar<'  discussed.  ! 


